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Abstract: The thrust distribution in e~^e~ annihilation is calculated exploiting its exponen- 
tiation property in the two-jet region t = I — T <^ 1. We present a general method (DGE) 
to calculate a large class of logarithmically enhanced terms, using the dispersive approach in 
renormalon calculus. Dressed Gluon Exponentiation is based on the fact that the exponenti- 
ation kernel is associated primarily with a single gluon emission, and therefore the exponent 
is naturally represented as an integral over the running coupling. Fixing the definition of A is 
enough to guarantee consistency with the exact exponent to next-to-leading logarithmic accu- 
racy. Renormalization scale dependence is avoided by keeping all the logs. Sub-leading logs, that 
are usually neglected, are factorially enhanced and are therefore important. Renormalization- 
group invariance as well as infrared renormalon divergence are recovered in the sum of all the 
logs. The logarithmically enhanced cross-section is evaluated by Borel summation. Renormalon 
ambiguity is then used to study power corrections in the peak region Qt ^ A, where the hier- 
archy between the renormalon closest to the origin (~ 1/Qt) and others (~ starts to 
break down. The exponentiated power-corrections can be described by a shape-function, as ad- 
vocated by Korchemsky and Sterman. Our calculation suggests that the even central moments 
of the shape-function are suppressed. Good fits are obtained yielding af^ (Mz) = O.llOitO.OOl, 
with a theoretical uncertainty of ~ 5%. 
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1 Introduction 



Event-shape distributions, as well as other observables that are not completely inclusive, 
tend to have significant contributions from soft and collinear gluon emission. Pertur- 
batively, this appears in the form of logarithmically enhanced terms, the well-known 
Sudakov logs. Another, more general source of enhancement of perturbative coefficients, 
which is associated with the running coupling, is infrared renormalons. These appear 
due to integration over particularly small momenta and therefore refiect the sensitive of 
the observable to large distance physics. The perturbative expansion of an event-shape 
distribution has quite a rich structure incorporating both sources of enhancement. Re- 
summation techniques exist for both Sudakov logs and renormalons, separately. Yet, in 
many cases, both are important. The Dressed Gluon Exponentiation (DGE) method 
advocated here takes into account both these aspects. 

It is understood how infrared renormalons are related to power-corrections through 
sensitivity to the phase-space boundaries [Q. It is also known that Sudakov logs appear 
due to the same type of sensitivity |0, ^ . However, the precise relation between the two 
phenomena was not addressed before, except in the specific case of Drell-Yan production 

3, where the emphasis was on the identification of the leading power-correction. DGE 
allows one to address this relation directly, and thus analyse the distribution in the peak 
region where the hierarchy between the renormalons starts to break down. 

Consider a generic event-shape variable y which is infrared and collinear safe , where 
the limit y — > corresponds to the extreme two-jet configuration. The logarithmically 
enhanced part of the differential cross-section has the form 0: 



where a = as{Q'^)/T^- When y <^ 1 a, naive fixed-order calculation is not adequate and 
the large logs must be resummed. However, this does not imply that the full perturbative 
calculation has to be carried out to all orders: the factorisation of soft and collinear radi- 
ation allows one to calculate the log-enhanced part of the cross-section, roughly speaking, 
"exponentiating the one-loop result" . 

It is useful to define the integrated cross-section 



so that -R(O) = 0. The logarithmically enhanced part of R{y) has the form R{y)\ 




(1) 




where the singularity in the two-jet limit, y 



0, is cancelled by virtual corrections. 



l + En=iK(l/)a" with 



2n 



Rn {y) = '^n,mL' 



m 
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where L = In ^ and rn,m = —Cn,m-i/fn- The full R{y) canQbe written as 

In R{y) = In R{y) \,^^ + In R{y) l^^^,^^ , (2) 

where the first term contains only logarithmic terms that diverge in the two-jet limit 
whereas the second is finite in this limit. 

Exponentiation means that In R{y)\^^^ has the following perturbative expansion, 

oo n+1 

lni?(?/)Log = E E Gn,^^L"^ a- (3) 

n=l m=l 

= Lgi(aL) + g2{aL) + ag^{aL) + o^g^^aL) + ■ ■ ■ 

where the function gi{aL) resums the leading-logs (LL) Gn,n+iL"'^^a"- to all orders, g2{aL) 
resums the next-to-leading logs (NLL) Gn-nL'^oP'-, etc. Exponentiation is a universal 
feature which follows from the factorisation property of QCD matrix elements (see e.g. 
[0), provided that the phase-space integration also factorizes. The latter condition makes 
the calculation of the exact coefficients Gn,m rather difficult. However, thanks to an 
intensive research effort in the last two decades, there is today a detailed understanding 
of Sudakov logs in a large class of QCD observables. In particular, in the case of event- 
shape observables in e^e~ annihilation, like the thrust, the jet mass parameters, the 
C-parameter, and jet broadening parameters, the state of the art is resummation of 
In I up to NLL accuracy [|, |, |]. 

Current phenomenology of event shape distributions (see |jlO| for a recent review) is 
based on NLL resummation, neglecting gk{aL) for A; > 3 in (^). The resummed cross- 
section is matched with next-to-leading order calculations that are available numerically 
ll|, [12|, [l^. This combination allows a good qualitative description of the observables. 



with the main caveat being the fact that "hadronization corrections", which are not 
under control theoretically, are required to bridge the gap between the calculation and 
the measurements. 

Before alluding to non-perturbative physics it is worthwhile to re-examine the per- 
turbative calculation. It was already shown in the example of the average thrust, that a 
major part of the discrepancy between the next-to-leading order calculation and the data 
can be explained by resummation of higher-order perturbative terms associated with in- 



frared renormalons |jT^. Physically the Single Dressed Gluon (SDG) resummation of |jT^ 
allows a better description of the branching of a single gluon emitted from the primary 
quark-anti-quark pair. Following the spirit of BLM |jl5|, 0, this resummation takes into 
account the fact that the physical scale of the emitted gluon is its virtuality, rather than 
the centre of mass energy Q, which is usually used as the scale of the coupling. As 



emphasized in [0, the essence of this approach is the idea that the perturbative expan- 
sion should be reorganized in a "skeleton expansion" in analogy with the Abelian theory. 
Then, the separation between running coupling effects and conformal effects is unique 
and physical. We shall see here, that in the specific context of the Sudakov region, one 



^There are various different ways to make this separation. The one we use here is motivated by the 
so-cahed log-R matching scheme. See [|| . 
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has access to sub-leading terms in the "skeleton expansion", corresponding to multiple 
gluon emission, through exponentiation. 

In order to examine the approximations in the standard NLL perturbative treatment 
of the thrust distribution, let us first recall two main features of the resummation (0). 
The first is the fact that the calculation of the functions gk{aL) is based on an integral 
over the running coupling. The second is that each gk{aL) resums a convergent series in 
aL (for aL -C 1). At first sight these facts seem contradictory: usually integrals over the 
running coupling translate into renormalon factorial divergence. The resolution is that 
the expansion (0) itself is divergent. The functions gk{aL) increase factorially with k, thus 
endangering the validity of the truncation at a certain logarithmic accuracy, e.g. NLL. 
Consequently, when examining the logs Gn,mL"^ at a fixed order n in eq. (H), one finds 
that the sub- leading logs are enhanced by a relative factor of ml compared to the leading 
ones. In the region where the perturbative treatment holds {aL ^ 1) this numerical 
factor can over-compensate the power of the log and thus invalidate the expansion. 

Another (related) argument for keeping track of sub-leading logs is the fact that in- 
tegrals over the running coupling yield at once all powers of L. Therefore, truncation at 
a certain logarithmic accuracy implies renormalization scheme and scale dependence. A 
renormalization-group invariant calculation of In R{y)\^^^ requires keeping all the logs. Of 
course, the exact calculation of all the logs is far beyond reach. However, a specific class of 
logs is fully accessible: these are the logs that are leading in /3o, G 

n,m ~ Po ^ 1 namely the 

ones associated with a single dressed gluon (SDG) emission. Within this class all the sub- 
leading logs, i.e. the functions gk{aL), are calculable. Moreover, since the exponentiation 
kernel is primarily associated with a single gluon, the entire logarithmically enhanced 
cross-section, which includes any number of gluons, can be generated in some approx- 
imation by exponentiating the SDG cross-section. The result is fully consistent with 
NLL accuracy, provided an adequate choice of the coupling (the "gluon bremsstrahlung" 
scheme [17|) is made. Our exponentiation procedure is based on the standard assumption 
that gluons emitted from the primary quark and anti-quark are completely independent 
and contribute additively to the thrust (this assumption holds for soft gluons). Given 
this assumption the SDG cross-section exponentiates under the Laplace transform at 
any logarithmic accuracy. Thus, through integrals over the running coupling we are able 
to calculate In R{y)\^^^ with the same formal logarithmic accuracy (NLL), yet avoiding 
any truncation of sub-leading logs, which would inevitably violate renormalization group 
invariance. In this way we resum at once Sudakov logs and renormalons. 

Comparison of the perturbatively calculated event-shape distribution with data re- 
quires some phenomenological model to deal with hadronization. Traditionally, hadroniza- 
tion effects are included using detailed Monte-Carlo hadronization models such as HER- 
WIG [|18[ or PYTHIA ||19||. An attractive alternative that became popular recently is to 
parametrise the effects in terms of a small number of non-perturbative parameters which 
control power-suppressed contributions. The basic assumption, which we adopt, is that 
the form of the most important non-perturbative corrections can be deduced from the 
ambiguities of the perturbative result |2T|, ^ ^, g. It was shown g, ^ § that the 
primary non-perturbative effect is a shift of the perturbatively calculated distribution to 
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larger values of y: 

1 da 
a dy 



. , 1 da 
a dy 



{y - AiA/Q) (4) 



where Ai is a non-perturbative parameter. According to the model of ||^, this parameter is 
expressed as an observable dependent factor which is perturbatively calculable times the 
"effective infrared coupling" (see The shift @) allows a fair description of the data 

at all measured energies and in a wide range of y values for several event-shapes in terms of 
a single parameter - the "effective infrared coupling" . The success of this perturbatively 
motivated approach suggests that the perturbative calculation itself, if pushed further, 
may be much closer to the data. It should be emphasised though that non-perturbative 
effects are essential in the two-jet limit when aL ~ 0{1), and in particular in determining 
the precise location of the distribution peak and its shape. In this region the above- 
mentioned approach, which uses a single non-perturbative parameter, is insufficient. In H, 



2^ , p8| , p9| , |30|| a more general approach was suggested, namely to replace the infinite sum 
of power-corrections of the form XnA"'/{yQ)^ which exponentiate by a single shape- 
function whose central moments are A„. 

As a consequence of infrared renormalons the resummation of In R{y)\^^^ contains 
power-suppressed ambiguities which signal the necessity of non-perturbative corrections. 
The fact that renormalons occur in the exponent strongly suggest that non-perturbative 
power corrections also exponentiate, in accordance with the shift and the shape-function 
approaches. We show that the renormalon ambiguity contains valuable information on 
the non-perturbative corrections, which is not restricted to the level of the leading power 
correction. In particular, it is possible to deduce the y dependence^] of each power cor- 
rection l/Q" from the corresponding residue of the Borel transform of the exponent. 
This allows us to make definite predictions concerning the form of the non-perturbative 
corrections. In the specific case of the thrust distribution the main conclusion from this 
analysis is that the even central moments of the shape-function are suppressed. More- 
over, in the case of the second central moment, an additional cancellation occurs between 
large-angle contributions, which are included in the shape-function of refs. ^|, and 
a collinear contribution which was not addressed before. 

Another important aspect which our approach incorporates is the complementarity 
between perturbative and non-perturbative contributions, namely renormalon resumma- 
tion and explicit parametrisation of power corrections |T^. It was clearly shown in |T^ 



that a consistent treatment of non-perturbative power-corrections as well as a reliable 
determination of requires renormalon resummation. Such resummation is performed 
here for the first time at the level of the differential cross-section. 

This paper deals with a specific observable, the thrust distribution. Nevertheless, 
DGE is a general method that can be applied in any case where Sudakov logs appear. 
It makes a direct link between Sudakov resummation and renormalons. This paper is 



^The idea that the dependence of power-corrections on external variables can be deduced from the 
functional form of the renormalon ambiguity is not new. It was used e.g. in ||2^ for the case of higher- 
twist in deep-inelastic structure functions and in |]3l| for the case of power corrections to fragmentation 
functions in e+e" annihilation. In the latter reference (see also [Q) it was justified by the 'ultraviolet 
dominance' assumption. Here we use this idea at the level of the exponent. 
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organised as follows: we begin in section 2 by recalling the calculation of the thrust dis- 



tribution in the single dressed gluon approximation ||3^. In principle, this calculation 
does not require any specific kinematic approximations, however, it is technically sim- 
pler if one neglects |]T^ 0, |3^ non-inclusive decay [^, |31| of the gluon into opposite 
hemispheres. Here we adopt this approximation and provide further evidence for its va- 
lidity at small values of 1 — thrust. We then observe that log-enhanced terms dominate 
the single dressed gluon result in a large range of thrust values. This justifies concen- 
trating on the log-enhanced cross-section of the form (|I|). In section 3 we calculate the 
resummed cross-section by exponentiating the logarithmically enhanced part of the SDG 
result. We also analyse there the renormalon structure of the exponent and the relation 
with power-correction models mentioned above. In section 4 we compare our approach 
with known results. First we refer to the NLL calculation which is based on solving an 
evolution equation for the jet mass distribution B. We show that the two calculations 



coincide up to NLL accuracy, provided an appropriate choice of the coupling is made |17 



Next we briefly recall the matching procedure of the resummed result with the full NLO 
calculation. In section 5 we study the phenomenological implication of the suggested 
resummation by fitting to data in a large range of energies and in section 6 we summarise 
our conclusions. 



2 Thrust distribution from Single Dressed Gluon 

The purpose of this section is to calculate the thrust distribution in case of a single 
gluon emission. The result will be used in section 3 to calculate the physical multi- 
gluon cross-section, through exponentiation. We begin by recalling the calculation of 
the Single Dressed Gluon characteristic function (sec. 2.1). Then, in sec. 2.2, we isolate 
the logarithmically enhanced terms and study their structure and finally, in sec 2.3, we 
construct the Borel representation of the result. 



2.1 Single Dressed Gluon (SDG) characteristic function 

In renormalon calculations [^, |2^, |I| the all-order perturbative result, which resums 
all the terms that are leading in /5o, can be expressed as a time-like momentum integral 
of an observable dependent Single Dressed Gluon (SDG) characteristic function, times 
an observable independent effective-charge: 

/ -^(e) Aeff(eg2) = / - [^(6) - ^(0)] -p{eQ') (5) 
JO e Jo e 

where Aik"^) = f3Qa{k'^) = /9otts(/i;^)/7r, and the bar stands for a specific renormalizat ion- 
scheme. In the large /?o limit (large Nf limit) A^k"^) is related to the MS coupling 

A{k^) = Poof {P)/7r by 

Going beyond this limit is discussed in sec. 4, where we identify A with the "gluon 
bremsstrahlung" coupling [jl^, thus making the DGE result exact to NLL accuracy. 
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The two integrals in (^) are related by integration by parts: 



and the function p(/i^) is identified as the time-like discontinuity of the coupling, 

The "time-like coupling" A^sifJ''^) in (^) obeys 

2 dAesifx'^) 



A ( -/i^ + ie) - A[ -/i^ - ie 



(7) 



(8) 



(9) 



For example, in the one-loop case A{k'^) = l/(lnA;^/A^) and 



11 /I fj.^' 

arctan — m 

2 TT U A2, 



(10) 



Eq. (^) which represents the all-order perturbative sum is not yet defined to power 
accuracy. The precise definition will be given in sec. 3.3 using Borel summation. A 
delicate issue, which was addressed in ^ |T^ is the fact that the dispersive 

integral (|^) over the time-like perturbative coupling (e.g. ([lOD) is formally well-defined 
but it differs from the Borel sum by power corrections, which are not entirely related 
to large distances. In the dispersive integral is given non-perturbative meaning by 
replacing the perturbative coupling by a non-perturbative coupling which is regular in 
the space-like region (i.e. it does not have a Landau singularity) and obeys a dispersion 
relation. Here we do not make this replacement. We refrain from attaching any non- 
perturbative meaning to (^. Large distance power-corrections will be deduced directly 
from the ambiguity of the Borel sum (sec. 3.5). Eventually, the normalisation of these 
power corrections can be expressed as an integral over an infrared-finite coupling, which 
is conjectured to be universal. 

The characteristic function J^{e) itself can be calculated either using the Abelian large 
Nf limit followed by naive non-Abelianization ||l|], or in the dispersive approach ||25|, |26|, 
2^ , referring to /i^ = eQ^ as the squared gluon mass. In case of inclusive observables, 
such as the total cross-section in e+e~ annihilation, the two calculations coincide. The 
two differ ||33|, |l^, ^ |30|, however, in the case of observables such as event- 



shape variables, that are not completely inclusive with respect to the decay products 
of the gluon. Here we shall use the massive-gluon characteristic function for the thrust 
distribution, which was calculated in p^. It corresponds to an inclusive integration 



over the decay products of the gluon, so that the possibility that a gluon decays into 
partons that end up in opposite hemispheres is not taken into account. It is important to 
emphasise, though, that the steps that follow can be repeated in very much the same way 
starting with the large Nf result, once available. In the case of the thrust, the inclusive 
approximation is expected to be good (some evidence is given below). Note that for some 
event-shapes, e.g. the heavy jet mass, this is probably not so. 
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Let us now briefly recall the calculation of the thrust distribution characteristic func- 
tion of ref. [0. At the SDG level, the distribution of t = 1 — thrust is 



1 da 

a dt 



{Q\t) 



C fi fif _ C r'^ df 

-f -A,s{eQ'me,t) = -f -p{eQ')[J^{e,t)-mt)]. (11) 
Po Jo e po Jo e 



The characteristic function JF(e, t) is obtained from the following integral over phase- 
space, 

jF(e, t)= / dxi(ix2 A^(xi, X2, e) (5 (1 — T(xi, X2, e) — t) (12) 

^phase space 

where CpciAi is the squared tree level matrix element for the production of a quark- 
anti-quark pair and a gluon of virtuality /i^ = eQ^, and 

' {xi + e)^ + {X2 + e e 



M{xi,X2,e) 



1 - a;i)(l - X2) 



;i-^i) 



X2 



(13) 



The integration variables Xi^2 represent the energy fraction of each of the quarks in the 
centre-of-mass frame. The energy fraction of the gluon is X3 = 2 — xi — a;2. 



We use the following |2^, ^ definition of the thrust 

y ^ Ei IPi-nrl _ Ei IPi-nrl 



(14) 



J:^E^ Q 

In case of three partons (a quark, an anti-quark and a "massive" gluon) it yields [|14 

1 — T(xi, X2, e) = min |l — xi, 1 — X2 , 1 — ^ {2 — xi — ^2)^ — 4e| . (15) 



Evaluating (|I2]) one obtains |^2[ the characteristic function for the thrust distribution 

^^(e,t)+ J-G(e,t) 



(16) 



where the dominant contribution J-'Q{e,t) corresponds to the phase-space regions where 
one of the primary quarks carries the largest momentum (T = Xi^2) and J^ci^i^) corre- 
sponds to the region where the gluon momentum is the largest. The superscripts / and 
h on J-'Q{e,t) denote low and high t values, respectively. These functions are given by 



^J(e,t) 



1 r 
7 



+ 



(l-t + e)2 + (l + e)2 
1 



In 



q-t 



2'- 

t 



t 



1 r 
1 



'l-t + e? + (l + ef 



In 



t{q-t) 



+ 1 



q 1 1 
1^2! 



q 1 1 
2- + -- 
t 2 t 



t-q 



+ 



^G(e,t) 
where q = 



' 1 1 
3- + 

, t q-t 



\-t 



^ + (l + e-g)2)ln^ + (2t-g) ("g + 1 + 



(17) 



;i-t)2 + 4e. 



2.2 Logarithmically Enhanced Terms 

Next we examine the SDG perturbative expression ([TT]) with the characteristic function 
of eqs. (|16D and (|^). For this purpose it is convenient to expand the time-hke couphng 
Acs{eQ'^) in terms of a fixed-scale space-hke couphng A^Q"^), 



m=0 i=0 ^ 



where, in general, K^j depend on the coefficients of the (3 function which controls the 
evolution of A^Q"^). Performing the integral ([TI|) over J-'{e,t) = — e^jF(e, t) order by 
order one gets the sum 

^E(^W')) i:i<m,hl{i) (19) 



SDG '^0 m=0 j=0 



where the functions hj{t) are the log- moments of the characteristic function, 

h,it)= / - In- ^(e,t). (20) 



e V e. 

The analytic calculation of hj{t) with the full characteristic function of eqs. (|T^) and (|T^ 
is rather involved. Alternatively, concentrating in the small t region, we distinguish in 
hj{t) between logarithmically enhanced terms which are not suppressed by any power of 
t and the rest, which we eventually neglect. 

As we shall see, the logarithmically enhanced terms, which have the form of eq. (0), 
can be easily computed analytically to all orders. Using (|T6|) for t < 1/3 we obtain 

"M-Ci{^nl)'M.ti + (.„1)VS(M) + (lni)Vj,(M). (21) 

We notice that a term of the form e"/ in JF(e, t) contributes to hj (t) a term proportional 
to (times a polynomial in In i) from the coUinear gluon emission limit of phase-space| 
where e = t, and a term proportional to t^n-fc (^^jj^ies a polynomial in In j) from the large- 
angle soft gluon emission limit of phase-space where e = t^. Taking the double expansion 
of (|T7|) at small t and e we have 

^S(M) - ^ + ^ (22) 

where we ignored higher powers of e multiplying a given power of 1/t, which are not rele- 
vant. From ( ^Tl) and (p2|) we conclude that only J-'q^g, t) contributes to the logarithmically 
enhanced cross-section and thus 

§A detailed analysis of the phase-space can be found in js^]. See figure 1 there. 



(23) 
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Performing the integral over e yields 



2 ^^^VT^^'^' - (j + l,L)- 2-^^' e^^r (j + 1, 2L) 
2- -L^+' -T -n 1 + 2"-'-') L' 



(24) 



where L = In j . 

Having evaluated the log-moments we can substitute them into (|1^) to obtain the 
logarithmically enhanced SDG differential cross-section to arbitrarily high order. The 
result (^) has a few features that deserve attention. First we note that the leading-logs 
in hj{t) depend on both collinear and large-angle soft gluon limits of phase-space, while 
the sub-leading logs are associated with the collinear limit alone. Next, we observe that 
the leading-logs in hj{t) originate in the e independent term 2/t in ^Q(e, t), while the 
sub-leading logs originate in the two e dependent terms, —e^/t^ — e/t^. This e dependence 
is reflected in factorial increase of the coefficients (see below) which is associated with 
the sensitivity of the thrust to collinear gluon emission. 

It is important to notice that the factorial growth is associated with sub-leading logs. 
Considering the different logs V' in a given log-moment hj{t) we find that with decreasing 
power k the numerical coefficient increases a,s j\/k\. Note, in particular, that the next-to- 
leading logs U (which are the only sub-leading logs taken into account in the standard 
analysis) have a small coefficient (|), while further sub-leading logs (which are usually 
ignored) are multiplied by large numbers. It is useful to analyse the square brackets in 
(^^ by considering two extreme limits for the function F{L) = e^T{j + 1,L). In the 
limit L — > 0, one finds factorial growth F{L) ~ j!. In the limit L — > oo one can keep 
only the leading logarithmic term in F{L), i.e. F{L) ~ L\ recovering the standard NLL 
result. The fact that the latter formally holds may be of little relevance in practice given 
that the entire perturbative treatment is adequate only when L <C 1/A{Q'^). One can 
now examine the validity of these two extreme approximations for phenomenologically 
interesting values of L. Consider for example the case t = 0.05, L = ln(l/0.05) ~ 3 at 
order j = 3: the value is F{L) = 77.8 while the NLL gives = 26.9 and the naive 
factorial gives j! = 6. Clearly, none of these extreme approximations is valid. One must 
retain all the logs. 

One may wonder to what extent it is justified to consider only the logarithmically 
enhanced part of the log moments. Factorial increase is certainly not restricted to the 
logarithmically enhanced part. The immediate answer is, of course, that the neglected 
part is parametrically smaller, by a factor of t compared to the one kept. This answer is 
not complete: as we know, parametric suppression can be compensated by large numerical 
coefficients. To be convinced that this approximation is valid at physically interesting 
values of t let us now compare the logarithmic part of the log-moments hj{t) with the 
full function which we evaluate numerically. Writing hj{t) = hj{t)\^^^ + Ahj{t), where 
the logarithmically enhanced part hj{t)\^^^ is given by ( p4|) and Ahj{t) are suppressed by 
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Figure 1: The functions hj{t), j = 1 through 4, in different approximations normahsed to 
the leading logarithm (LL). The complete result (calculated numerically) is the full line 
and the sum of the logarithmically enhanced terms hj{t)\y^^ (calculated analytically) is 
dashed, the NLL approximation is dotted. For the first moment, the jSo dependent part 
of the NLO result is also shown for comparison. 



a relative factor of t, the first few log- moments are 
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The hjit) and thus Ah jit) are calculated by integrating (^) numerically where jF(e, t) 
is the derivative of the full characteristic function (0). The comparison for the first four 
log-moments is presented in fig. |1|. We see that while the NLL approximation is quite 
far from the actual values of hjit) even at rather small t, there is a fairly wide range in 
which the logarithmic approximation hjit)\^^^ is good. At large t values the gap between 
the two increases. However, since our interest here is in improving the calculation in the 
small t region, we consider this agreement satisfactory. It should be kept in mind, though, 
that our resummation procedure is not the appropriate way to improve the calculation 
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of the cross-section at t ^ 0.2. We comment that for higher log-moments j > 5 the 
approximation deteriorates even at smaller t values; however, the resulting discrepancy 
at the level of the SDG cross-section ([19|) is not large. We shall revisit the significance 
of terms that are suppressed by t in the phenomenological analysis in sec. 5. 

Finally, using (j^), the first few orders of the SDG calculation (19) of the log-enhanced 
cross-section are given in the MS scheme by 



1 da 
a dt 



SDG 



{(2 L - 1.5) A(Q^) + (3 + 1.8333 L - 3.75) A{Q^f (25) 

+ (4.6667L=^ + 8.5000 - 8.5242 L - 5.6481) A{Q^f 

+ (7.5000 + 21.833 - 15.859 - 40.585 L + 2.025) A{Q^)^ 

+ (12.400 + 48.500 - 29.339 - 175.72 - 47.164 L + 40.22) A{Q^f + ■ • ■} 

The leading and next-to- leading logs are consistent with previous calculations 0. 
In addition, at NLO we can compare the next-to- next-to-leading log (—3.75 in ([25|)) 



with the exact numerical computation |jT3[. To make an explicit comparison, one should 
separate the contributions of different colour factors in the NLO coefficient (this was done 
before using the EVENT program [jl2[ in ||4^, ^]) and then fit the numerical results to 



a parametric form which includes a log-enhanced part, and terms which are suppressed 
by powers of t. The result of such a fit, where the first two logs are fixed to their exact 
values, gives the following log-enhanced NLO expression in the large j3o limit. 



a dt 



{(2 L - 1.5) A(g2) + (3 + 1.8333 L - 3.77 ± 0.07) A{Q^f} .(26) 



Details concerning the fit and the error estimate can be found in the Appendix. The 
closeness of the NNLL coefficients is encouraging. It strongly suggests that the impact 
of the non-inclusive splitting of the gluon into partons that end up in opposite hemi- 
spheres |HT|, which is not taken into account in our calculation, is small, at least at 
the logarithmic level. A similar comparison, which goes beyond the logarithmic level, has 
been performed for the average thrust in |T^, ^ and the higher moments (t") in p2[ . 
It was found that the non-inclusive effect is rather small for the average thrust (4.4%) 
while it increases for the higher moments. This implies that the effect is mainly impor- 
tant in the large t region. Being interested here in the small t region and in particular 
in the logarithmically enhanced part of the cross section, we shall neglect this effect al- 
together. The reader should keep in mind, however, that the evidence provided above 
concerns directly only the Abelian part of the coefficients. Only through the assumption 
that running coupling effects dominate (making the large jSo approximation legitimate) 
this evidence becomes relevant for the no n- Abelian non-inclusive effect. Moreover, in the 
context of power-corrections within the infrared-finite coupling approach, the impact of 
the non-inclusive effect in the non- Abelian case was shown to be more significant than in 
the Abelian one [BUI, HD . 
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2.3 Borel representation of the SDG result 



A convenient way to deal with perturbation theory to all orders is Borel summation. In 
this section we construct the scheme-invariant Borel representation of the logarithmically 
enhanced part of the SDG result. This will be used in the next section to calculate the 
multiple emission cross-section. We start with the scheme-invariant Borel representation 
of the coupling |^ ^ 



A{e) = dz exp (-zlnitVA^) Ab{z). 



(27) 



In the case of one-loop running coupling (2;) = exp zlnA^/A^^. Taking the time-like 
discontinuity we get 



Aff(Q') = dz exp (-zlngVA 



2^ Ab{z). 

TCZ 



(28) 



In eq. (0) we use the part of jF(e, t) that contributes to log-enhanced terms (as in (p3D), 



a at 



* de ^ 



e 



e 



(29) 



Substituting Ac{[{eQ'^) from (^) and integrating over e we obtain the Borel representation 
of the logarithmically enhanced SDG cross-section, 



1 da 
a dt 



{Q\t) 



/?0 



1 / 

-/ dzBsj,G{z,t) e^^[-z\YiQ'^/K' 
t Jo ^ 



sm Tiz 



HZ 



Aniz) 



with 



Bsvg{z, t) = - exp {2zL) 

z 



1-z 2-z 



exp (zL) 



(30) 



(31) 



where L 



Ini 



Note that the first exponent in the Borel transform is associated with 



the large-angle soft gluon limit of phase-space, while the second is associated with the 
coUinear limit. The Borel integral is well-defined at the perturbative level thanks to the 
cancellation of the pole aX z = between these two contributions. Moreover, provided 
that the coupling does not enhance the singularity (e.g. the case of one- loop coupling), 
the poles at ^ = 1 and z = 2 are regulated by the ^^^7^ factor, which originates in the 
analytic continuation to the time-like region. 

Note that the Borel integral converges at infinity only provided that 2L < InQ^/A^, 
or > A, reflecting large-angle gluon emission sensitivity. A weaker condition appears 
from the collinear limit: L < InQ^/A^, or Q'^t > A^. These convergence conditions 
signal the fact that the small t region is not under control in perturbation theory, and 
that power corrections of the form A"/ (Qt)" (and eventually also A^"/ (Q^t)") should be 
included [E^, 120, pO]. We will return to this point in more detail in the next section. 
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3 Exponentiation 



The calculation of the previous section concerns a single gluon emission from the primary 
quark-anti-quark pair. Close to the two-jet limit, resummation of multiple soft and 
collinear gluons is essential. This resummation is achieved by exponentiating the SDG 
cross-section. As mentioned in the introduction, our basic assumption is that emission of 
several soft or collinear gluons can be regarded as independent. The total contribution to 
the thrust is then computed as the sum of the individual contributions of these emissions: 
t = J2k=i ^k- For large- angle soft gluons, which give the dominant contribution at small t, 
additivity holds since the contribution of each gluon to t is proportional to its transverse 
momentum. This assumption fails, of course, for sufficiently hard gluons. In general, the 
contribution of several gluons to the thrust, as any other event-shape, are correlated. Such 
a correlation emerges from the modification of the phase-space limits in the integration 
over the gluon momenta (the quarks do recoil) as well as from additional interactions 
between the emitted gluons. These correlations are expected to be important only at 
large values of 1 — thrust, away from the two-jet configuration. Since they may have an 
effect at the logarithmic level, our calculation of sub-leading logs (beyond NLL accuracy) 
that are sub-leading in /3o[|| is only partial. It remains for future work to quantify the 
effect of these correlations. Here we simply assume that it is small. 



3.1 Exponentiation Formula 



The result of the previous section, ^^{Q'^,t) , has a probabilistic interpretation as 
the probability distribution, of the random variable t G (0, 1), given that only one gluon 
is emitted. Since, in reality any number of gluons n can be emitted, we sum over these 
possibilities 



1 da 
a dt 



{Q\t) 



oo n 
n=0 "'■ k=l-^^ 



dtk —-jri^k] 
a at 



6 t 




m 



(32) 



where the singularity due to zero gluon emission probability = is cancelled by virtual 
corrections, S{t). The combinatorial n\ is introduced to avoid multiple counting - the 
expression is symmetric for any permutation of the n gluons. Using the Laplace conjugate 



variable z/, the S function can be written in a factorized form as (see e.g. |22 



" f dv 

^(t-E4) = y^^exp 




(33) 



where the contour C runs parallel to the imaginary axis (to the right of all the singularities 
of the integrand, if any). Changing the order of the v integration and the sum over n in 

^Such correlations do not affect the logs that are leading in /3o. These terms are not modified by the 
exponentiation and their calculation in the previous section is exact barring the inclusive approximation 
made in the definition of the thrust. 
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(P^ we get the exponentiated form, 



1 da 
a dt 



dv ^ 1 \ Ida 



c Ztti jo a dt 



a dt 

{Q\i) 



- 1 ) dt 



-vt 



l]di 



Defining 



we have 



5pT((5^,lnz/) = (e 



1 da 



-vt 



{Q\t) 



a dt 

and finally, integrating over t, 

R{Q\t) 



1 1 da 



(e'"'* - l) (it, 



SDG 



c 27ri 



e exp 



dv 

c 2TTi V 



e''* exp 



Note that the terms which are leading in /5o in the coefficients of j^) 

log-enhanced part of the exact result, are the same as the ones in 
coefficients correspond to a single emission. 



(34) 

(35) 
(36) 

(37) 
or in the 



These 



SDG 



3.2 The Exponent in the Conjugate Variable 

Starting with the Borel representation (^) with the Borel function (^) we now calculate 
the exponent as defined in (|35D, 

= dzBJz)exp(-z\nQyA 
>G do Jo ' \ 



9 , . n\ Sm TTZ -r , , 



(3J 



with 



Bjz) 



1 dt 
T 



z 1 — z 2 — z 



(39) 



To evaluate B^{z) we use the integral 

JO t ^ ^ z 



(40) 
(41) 



where 

7(-^, z/) = r(-2;) - r(-2;, z/). 

Note that (^OD is regular at z = due to cancellation between the 7 function and the 
simple pole - the first originates in the e"'^* part, i.e. in real contributions, while the 
latter in the unity, namely in virtual corrections. We thus find 



Bjz) 



e'^'-'^^{-2z,u) + - 



2 1 1 
- + + 



z 1 — z 2 — z 
14 



e^^-'^^i-z,u) + - 



(42) 



Similarly to (piD, this Borel function is well defined at the perturbative level thanks to 
cancellation of the pole at 2; = between the two terms in (^21); the first corresponds 
to the large-angle soft gluon emission and the second to collinear gluon emission. In 
contrast[|] with BsoGi^^i) in (0)) B^{z) has a rich renormalon structure, making the 
integral (^) ambiguous at the level of power accuracy. Note that 7(— -z, z^) has the same 
singularity structure as r{—z): a simple pole at any non-negative integer. It follows that 
By{z) in (]4^ ) has two sets of infrared renormalons: the first at all positive integers and 
half-integers due to large-angle emission sensitivity and the second at all positive integers 
- due to collinear sensitivity. In the latter, the renormalons at 2; = 1 and z = 2 become 
double poles. It is important to note that the sin(7r2;) factor in (^) regulates simple 
poles and transforms double poles into simple poles at all integer z values, while it does 
not affect the singularities at half integer z values. We will return to this point in the 
context of power-corrections in sec. 3.5. 

An alternative representation of (e~'^*)gpQ in terms of a characteristic function can be 
obtained using eq. (P^), namely 



-ut\ 



Jo Jf^ e 



Po Jo e 



e 



e 2 



(43) 



where the order of integrations over t and e has been changed. The characteristic function 
is given by. 



2 + z/e 



1 2 2 
- z/ e 



+ 



+ 



2 2 



Ei(l, v^z/) + 



-V e 



2+ -z/^e^ 



Eid, z/e) 



+ 



21n(v^z/) -21n(z/e) + - e + - -. 



3 1 

-z/e 

2 2 

3 



(44) 



Given a modelQ for the coupling which is regular in the infrared 



I , the integral ( ^3]) 

could be readily evaluated. We prefer, however, to complete the perturbative treatment 
without any additional assumptions on non-perturbative physics. For this purpose we find 



the Borel representation (|42|) most suitable. In the next section we calculate this integral 
explicitly, but before doing so it is worthwhile to examine more closely the structure of 
the sum, term by term. 

Since we are interested in large v we can replace 7(— 2;, z/) in (42) by r(— 2;). The 
logarithmically enhanced terms are not modified by such a replacement. Expanding the 

"At first sight it seems surprising that the Laplace-transformed By{z) has a completely different 
renormalon structure from _BsDG(2^ii)- The reason is, of course, that B^{z), thanks to the integration 
over t, is sensitive to the small t region, where the exponents exp {2zL) and exp (zL) in are large 
and impose convergence constraints on (|30|). 

*An example is provided by B^. 
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integrand in 



in powers of z, we obtain 
(In z/)" 



E 

n=0 



2 (i-nc„ 

m>max{l,n— 1} 



^m— 1 



(45) 



m>max{n+/,n+l} 



m— 1 



Z>0 



where the first term in the square brackets originates in the 2/t part in J-'{e,t) from 
both the soft and colhnear hmits of phase-space, while the second originates in the e 
dependent terms in JF(e, t) from the colhnear limit alone. The numbers are defined 
by T{~z) = —J2T=-i CkZ^. This means that c_i = 1, Cq = 7^, Ci = 7r^/12 + 7|;/2, etc. 
Higher contain higher Q numbers {i <k + 1), yet numerically (A; > 1) are all close 
to 1. 

Next, performing the Borel integral term by term we get a series in the coupling A{Q'^) . 
It is then straightforward to rewrite the sum as 



SDG 



^E^wt"'/.(^(g^)in^) 

fc=i 



(46) 



and perform the summation over all powers of ^ = A{Q'^)\nv. The functions fk are 
given by 



{fc-l)/2 



U0= E 



-vr 



2\M 



.=0 (2^ + 1)' 



where 



fc-2 



Mcol(fc) 

Mia(A;) 



2cfc_2+ E(l + 2 



— m— 1^ 



(47) 



(4J 



m=0 
Cfc_2 



and 



= r ^^"^"^ Hypergeom([l, - 2 + no], [1 + no], 

r(l + no) 



with no = max(0, 2-u + 3 — /c). In (|48|), Mcoi(/i;) and Mia(/i;) correspond to the colhnear 
and large-angle contributions to the Borel function, respectively. Explicitly, this yields 



/i(0 = 2(l-01n(l-0-(l-201n(l-20 
/2(0 = -2 7i.(ln(l-0-ln(l-20)- ^ln(l-0 



(49) 



and 



/3(0 = 0.804/(1 - - 2.32/(1 - 20 UO = 0.779/(1 - 0' - 2.68/(1 - 20' 

/5(0 = 2.32/(1 - 0' - 5.00/(1 - 20=^ h{0 = 6.12/(1 - 0' - 15.32/(1 - 20' (50) 

/7(0 = 23.8/(1 - 0' - 61.12/(1 - 20^ /8(0 = 120.9/(1 - 0' - 305.52/(1 - 20^ 
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etc. Here the function /i resums the leading- logs (LL), /2 resums the next-to-leading logs 
(NLL) and so on. The effect of the analytic continuation to the time-like region is taken 
into account through the summation over u in (0). Note that we omitted a constant 
term in eq. (|50|). Such a constant does not necessarily exponentiate and it can be shifted 
into the "remainder function" R{y)\noa_iog (0)- 

The first two functions /2(0 already known from When comparing 

/2(^) in (|49D to one should be careful concerning the definition of the coupling. Strictly 
speaking, our calculation is done in the large /Sq limit, so jSi terms as well as other terms[[] 
which are sub-leading in (3o are missing. However, these terms can be generated upon 
replacing the large /3o couphng, which is one- loop, by a two-loop coupling in the "gluon 
bremsstrahlung" scheme (see sec. 4) inside the integral (0). Such a replacement will 
be made in sec. 3.3 and then in the phenomenological analysis in sec. 5. 

The functions fkiO^ k > 2 are singular at ^ = 1/2 and at ^ = 1; the first is 
associated with large- angle emission uA/Q ~ 1 and the second with coUinear emission 
uA'^/Q'^ ~ 1. The singularity (at both limits) becomes stronger as the order k increases, 
in particular fk ~ 1/(1-2^)''"^ close to ^ = 1/2 and fk ~ 1/(1 -^''"^ close to ^ = 1, for 
any k > 3. Since soft gluon resummation has a universal structure, we expect that this 
singularity structure in ^ will be common to different physical processes. An example is 
provided by a recent calculation of NNLL for deep-inelastic scattering and the Drell-Yan 
process [Q. There, the singularity structure of fk{0 for /c = 1 through 3 is indeed similar 
to (|49|) and (^Ol). In particular fsi^^) has simple poles. 

As usual in perturbation theory the expansion (^) is divergent. In principle it can 
be given sense e.g. by truncation at the minimal term. Note that this truncation is 
cumbersome since the minimal term will correspond to different orders depending on the 
value of u. The standard procedure amounts to truncating this expansion at NLL, so 
fkiO k > 3 are just neglected. To see whether this is legitimate let us now examine 
the numerical values of fk{0- 

The functions fk{0 plotted in fig. 0. It is clear from the plot that the expan- 
sion (^) is divergent. In particular, the ratios /fc(0//fc+i(0 increase with k, reflecting 
factorial growth. The actual effect of these increasing coefficients depends on the value 
of the expansion parameter, here A^Q"^) = Poa^Q^). A{Q^) is defined in the "gluon 
bremsstrahlung" scheme [0 (the reason for this choice is explained in section 4). In 
fig. we show the individual terms in (|i6| ) at two representative centre-of-mass energy 
scales: Q = 91 and 22 GeV. First note that the NNLL (/s) gives a significant contri- 
bution in all cases. It should be emphasised that the leading part of its contribution is 
taken into account in the standard NLL calculation in the process of matching with the 
NLO result. Next note that N^LL (f^) also gives quite an important correction: ~ 8% at 
Mz and ~ 18% at 22 GeV. A closer examination shows that the hierarchy between the 
contributions of sub-leading logs depends strongly on both the energy (the value of the 
coupling) and l/u. At sufficiently large l/u values, i.e. far enough from the two-jet limit, 
the expansion (^61) appears to be convergent at the first few orders. Consider for example 
the case 1/z/ = 0.075 at 22 GeV: the contribution of (full line) is smaller than that of 

^Such terms originate [p| in the singular part of the NLO sphtting function. 
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Figure 2: fk{Cj for = 1 (LL) through A; = 8 (N^LL) as a function of ^ = A{Q'^) Inz/. 
The absolute value is plotted on a logarithmic scale. 

/a and /s (dotted) is yet smaller. But then the trend changes: /e (dashed) contributes 
slightly more than /s and the series starts diverging. In this particular example the con- 
tribution at the minimal term is about 7.8%. Thus a resummed perturbative calculation 
should not be expected to be better than that. For smaller 1/v values, and likewise at 
lower energies, the expansion diverges faster: the order at which the minimal term is 
reached becomes smaller while its contribution becomes larger. At some point, close to 
the distribution peak, the expansion becomes meaningless. Technically, this behaviour is 
a direct consequence of the enhanced singularity of the functions fk and their factorially 
increasing numerical coefficients. This is the way infrared renormalons show up in this 
framework. Physically it is well expected that the perturbative expansion would gradu- 
ally become less predictive and eventually irrelevant as a quantitative evaluation of the 
cross-section, since it misses non-perturbative soft gluon contributions or, in other words, 
hadronization effects. We shall return to this issue in sec. 3.5. Until then, our goal is to 
complete the perturbative evaluation of the cross section. To achieve this we leave the 
expansion ( ^6D and revert to Borel summation ( pSj) which allows us, in rough terms, to 
resum the renormalons first, and then complete the summation over the logs. 
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Figure 3: The relative significance (in absolute value) of terms A{Q'^)'^ ^ fki^A{Q'^) Inu 
in (ID for = 2 (NLL) through k = 7 (N^LL) as a function of 1/z/ (1/z/ ~ t in the 
small t region). In both plots, the terms are normalized by the LL function (1 on the 
vertical axis). 



3.3 Evaluation of the Borel sum 

Our next task is to evaluate explicitly the Borel sum (PS]). An exact analytic integration 
is difficult. However, it is possible to make a systematic approximation of the Borel 
function and then perform the integration analytically. 
We begin by writing the Borel function By{z) as 

B^{z) = -V{2z) - (2 + 1 + J_\ v(z) (51) 
z \z I — z 2 — zj 

with 

V{z) = e'^'"'T{-z) + (52) 

Since V{z) has only simple poles at positive integer z values, it is natural to approximate 
it by a rational function. Clearly, we would like the approximation to coincide with the 
function at least at the first few orders in the expansion around z = 0. One possibility is 
then to use Fade approximants based on some truncated sum. However, since we know 
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the singularity structure, it is better to use an approximation of the form|^ 



Pmax 



y(^) - E rfr, (53) 
p=i y ^ 



where are set such that the expansion of this approximant coincides with that of (|52D 
up to order 2;^max-i_ This means that r^, for every 1 < p < Pmax? is a polynomial of order 
Pmax in In i^, 

Pmax 

rp=i: Pp,, (In^)'^ (54) 

where Pp^q are (irrational) numerical coefficients. 

Next, B^{z) can be written as a sum of simple and double poles, 

Pmax Pmax ^Col 

i?.(-)-E4^+E^-7T^-7^ (55) 

- 2; ^1 p - 2; (1 - z) (2 - z) 

where f^^ and fp°^ are linear combinations of {rpY'^^, corresponding to the large-angle 
(first term) and coUinear (second term) parts of (|^, respectively. Simple poles appear 
at all positive integers and half integers 1/2 < p < Pmax/2 due to the large-angle part, 
and at all integers 3 < p < Pmax due to the collinear part. The double poles at 2; = 1 
and z = 2 appear due to the collinear part. As an example, consider the case of Pmax = 3 
where one has 



2ri ^ -?,ri- f rg ri - rg , |r3_ _ ri rs 

z V-z ^- z 2- z -6- z {l-zf (2 - 2)2 ■ 



2 ^ 1 ^ 2 



Finally, given a specific coupling and the function (|55D the Borel integral ( pSf ) can be 
evaluated analytically. The simple pole integrals are of the form 



/„(s) = / c/^expf-^lngVA^) ^ Ab{z) (56) 

JO ^ ' p — z nz 

= / a^; exp (— A.b\z) 

Jo p — Z TIZ 



where Ab{z) is the Borel transform of the MS coupling and s = InQ^/A^. Regarding p 
as a complex number, Ip{s) can be written in terms of Ip{s), 

-^p('^) = / dzexp(—zs) Ab(z), (57) 

JO p — z 

by analytically continuing in s, Ip{s) = ^ [Ipis — in) — Ip{s + ivr) j, and 



W = - 
p 



Ip{s) - lim Ip{s) 

p=0~ 



(58) 



similar approximation for the Borel function was used in in a different application. Comparison 
with Pade approximants showed consistent results. 
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The double pole integrals 



poo 

/ dz exp 
Jo 



-zs 



{p - zY 



TCZ 



Ab(z) 



(59) 



can be simply obtained using I^^\s) = ~-^Ip{s). 

In case of infrared renormalons in space-like observables Ip{s) should be evaluated for 
real s and p where it is ill-defined. Then Cauchy principal value can be defined taking 
the real part of ([57|). Our definition for the time-like renormalon integral is motivated 
by this regularisation. Here there is no need to take the real part: thanks to the analytic 
continuation procedure, ( [F7| ) is evaluated above and below the cut in a symmetric way 
(the relevant complex parameter is ps), yielding a real value for Ip{s) and I^\s). The 
ambiguity of these quantities is directly related to the ambiguity in Ip{s). The definition 
of the Borel-sum through the space-like renormalon integral also allows to introduce a 
cut-off regularisation . 

Ip{s) can be written in terms of known analytic functions. The simplest case is the 
one-loop running coupling, where Ab{z) = 1. Here we have 



Ipis] 



one— loop 



— Ei(l, —ps) e 



■ps 



(60) 



The more difficult two-loop case was dealt with in |[T^ (see sec. 3.4 there) using results 
from p2|. It is useful to express the coupling using the Lambert W function P5[], 



-loop 51 + w{s) 



(61) 



where 5 = /3i/ /3q and f\w{s) = W^i ^— e ^j. In these variables the two-loop renor- 
malon integral (^^ is simply 



two— loop 



- [w{s)p6Y^ e^^^'^P^ r {-p6, w{s)p6) . 



(62) 



Altogether the Borel integral (|38D is given by 



f3o 



Pmax 



n=l 



'^n/2-^n/2( 



Pmax 



-col 



I^(s)-nli'\s)-r2l'2'{s 



(1)/ 



p=i 



(63) 



where the (polynomial) dependence on Inz/ is fully contained in the coefficients fp and 
Tp, given by eqs. (|52|) through ([5^), while the dependence on the external energy scale 
s = InQ^/A^ is fully contained in the renormalon integrals Ip{s). Eventually, the log- 
resummation order Pmux should be set large enough, so that higher orders can be safely 
neglected. In sec. 5 we show that in practice Pma.x — 6 is enough to make the result stable. 



§The analytic continuation in the calculation of Ip{s) should be done carefully: one should use the 
correct branch n of the Lambert W function, Wn- In practice this is simple since for Nc = 3, Nf < 6 
the time-like axis is always contained in the n ~ ±1 branches (see |^ for more details). Also note that 
with the current definition of the two-loop coupling the Landau branch point is at — A^. 
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3.4 Evaluation of the inverse Laplace transform 



The last stage of the calculation is to evaluate the inverse Laplace transform (P7|), 

R{Q',t) = I -^e'^*exp{5p,(g^lnz/)} (64) 

res Jc Zm V ^ J 

with 5pt((5^, Inz/) = (e^'^*)gj3Q. We apply the same technique used in (see section 
5 there). However, at a difference with the latter, we must keep track of sub-leading 
logs so we refrain from making any approximation. Defining u = vt and expanding 
SiQ"^ , In I/) = S{Q^, L + Inw) (where L = In j) in powers of \nu, one has 

R{Q\t) = f -^exp\u + S{Q\L + lnu)} (65) 

res Jc 2m U ^ ^ 

/ exp n + rfl(Q^ L) In n + ^ tt4(Q', i:)(ln uf) 
Jc Zirt 11 I ^ 2 I 



where dk{Q'^,L) = -^S{Q'^ , L). Finally, we calculate R{Q'^,t)\^^ by expanding 

exp|f:^4(Q^^)(ln«)'=)| 

in powers of In u and using the integral 

/c2^^^"")"^"P^" + ^^""^ = |^f(T^- ^''^ 

Our final expression for R{Q'^,t)\^^^ includes several special functions: derivatives of 
the r function from (|66|) and the expressions for Jp(lnQ^/A^), written in terms of the 
exponential integral function in the one- loop case (|60D, or the incomplete F function and 
the Lambert W function in the two-loop case (^). We therefore chose to evaluate it 
using the computer algebra program Maple. 



3.5 Power- Correct ions 

Since hadronization is not described by perturbation theory it is a priori clear that there 
are non-perturbative corrections to the resummed cross- sect ion. In this section we study 
the structure of these corrections. We already saw that the perturbative sum is, in 
general, ambiguous due to infrared renormalons. Having used a specific regularization of 
the renormalon singularities we are bound to address the question of the related power- 
suppressed ambiguity. The appearance of infrared renormalons does not only signal 
the existence of non-perturbative power-corrections, but also suggests their functional 
dependence on the energy and the event-shape variable. One should not expect, of course, 
that the entire non-perturbative contribution would be deducible from a perturbative 
result. Appropriate non-perturbative parameters will be introduced and eventually fixed 
by the data. Upon including these non-perturbative corrections the result becomes well- 
defined to power accuracy. 
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If the Borel sum (^) is supplemented by a non-perturbative contribution, eq. 
becomes 

R{Q\t) = [ -^e^*expf5pT(Q^lnz/)} e^^ {S^^{Q\\nu)] . (67) 
Jc I'm V ^ > ^ > 

In sec. 3.4 we evaluated the inverse Laplace integral ignoring the non-perturbative part. 
Thanks to the factorized form of (^) non-perturbative corrections can be included by 
convoluting 0, 12^, ESI, EUl the perturbative result (K|) with a non-perturbative function: 



R{Q\t) = C dtR{Q\t-i) U^{Q\t) = R{Q\t) ^ UiQ^t) (6. 

JO I'^s res 



where 



U{Q',t) = J^^e-"' exp{s^,{Q\\nu)} . (69) 

In order to learn about the functional form of iSnp(Q^, In i/), we shall now examine 
the renormalon structure of iSpt(Q^, In i/). Considering the integrand of ( pSj ) with By{z) 
of eq. (1^) (see description following the equation), we identify potential singularities 
at all positive integers and half integers. On the other hand, the sinvr^; factor which 
originates in the analytic continuation to the time-like region introduces an explicit zero 
at all integer z values. It is important to realize that the complete cancellation of some 
of the singularities may be an artifact of the large /3o approximation. A similar situation 
occurs in any time-like observable, e.g. in the total cross-section There, having 
an operator product expansion, one can relate the singularity strength of the Borelp| 
pole 1/ {p — z)^^*^ to the anomalous dimension of the corresponding operator (s) Op of 
dimension 2p: k = pPi/Pq + 7o//?o, where [d/ d\n ^"^ + •y] Op = with 7 = 70 («s/vr) + . . .. 
Thus, in general, a large-/?o calculation does not yield the correct singularity strength: 
not only the residue, but also k is expected to be modified when sub-leading terms in 
/3o are included. Being aware of this fact (and having no operator product expansion to 
supply additional information) it is reasonable to allow for modifications of order 1 in 
the singularities of B^{z), 



z 



where we assume \Kp\ < 1. Thus, instead of simple and double poles (m = 1,2), the 
Borel integrand has cuts which survive the suppression by the analytic continuation and 
introduce some ambiguity into Sp^^Q^ ,hiu). Note that, we will not consider the possi- 
bility of a modified u dependence of the exact Borel function compared to the calculated 
one. We allow only a modification of the overall coefficient and the singularity strength. 
Consider first a scheme-invariant Borel integral in the space-like region, 

J(s) = / e-'%z)dz, (71) 



^Note that the pPi/ 0^ term appears in the standard Borel transform but not in the scheme-invariant 



one |36 
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where the Borel function behaves as (^) with m = 1 near z = p. The Borel integral can 
be written as a sum of two terms, 



lis) 



dz e 



(p-z) 



dz e 



l + K 



{-sY 7(-/€, -sp) e-'^ - {-sY r(-/€) e 



sp 



i-sYTi-K,-sp) e 



sp 



Next we use the fact that the ambiguity A/(s) is purely imaginary and originates uniquely 
in the second integral where z > p, i.e. in the complete T function part[J]. Taking l/vr 
times the imaginary part, we obtain 



Ai{s) 



sm VTK 



TT 



sT(-/t) e 



-sp 



T(1 + k) 



-sp 



(72) 



The ambiguity for the m = 2 case can be calculated using 

= -j-i{s) : one finds 

that A/(i)(s) = A/(s)s. Substituting s = \nQ^/A^ we see that the ambiguity A/(s) 
scales as a power (Q^/A^)^, up to a logarithmic factor, (inQ'^/K^^ , which remains 
undetermined as long as k is unknown. Next, the ambiguity of the corresponding Borel 
integral in the time-like region 



Us) 



_~ smiTZ 

e Hz) dz, 

nz 



can be deduced by analytically continuing (|7^) in 
taking times) the imaginary part. 



-Q^, i.e. s 



(73) 

s ± ZTT, and 



AJ(.) 



1 (g' + TT , 

Tip r(i + K 
1 1 



-sp 



TT 



COS K arctan — ) sin(7rp) — sin k arctan — I cos(7rp) 



TT 



Tip r(i + k) 



KTC 



sin(7rp) cos(7rp) 



(74) 



where the expression in the second line is obtained expanding in 1/s and neglecting 
terms and higher. Similarly, we obtain 

1 1 



A/«( 



e '^^ [s sin(7rp) — (1 + K-)n cos(7rp)] . 



(75) 



TCp r(l + K,) 

Given the general singularity structure of the SDG Borel function By{z) around z = p, 



b{z) = ^ + 



(76) 



our result for the ambiguity ((0) and (|75D ) suggests the following assignments of non- 
perturbative contributions. 



Tip r(l + 



r'ph + rp) sin (up) - ( (1 + /tp) r^^^ + Kp Vp/ s) n cos (np) 



(77) 



"a similar mathematical analysis appears in |]l4| in the context of a renormalon integral with a 
two-loop running coupling. 
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where s = InQ^/A^ and p2p are dimensionless phenomenological parameters to be fixed, 
eventually, by the data. On general grounds we expect that p2p ~ It is important 

to keep in mind that the actual values of p2p depend on the regularization prescription 
used for the Borel sum. For example, when the perturbative sum is defined using a cut-off 
(on the Euclidean momentum), p2p can be interpreted as moments of an "infrared finite 
effective coupling" Aj^p{k^) [T^, which is assumed to coincide with the perturbative 
coupling y4pT(fc^) above some scale pj (A <C /xj -C Q), 



P2p 




(78) 



where 



m;^p'j) 




(79) 



Here we define the perturbative sum a la Borel. In this case p2p is 



P2p 




^0) 



where 



m-(/x?) = P.V. 



Ml 



P 



1^ 




and P.V. indicates the principal value prescription. Note that in this case p2p are in- 
dependent of Pi. a detailed comparison of the two regularization prescriptions can be 
found in [|l^. The explicit expressions for mp'^(/i|) in the one-loop and two-loop cases 
can be obtained by substituting s — > ln/i|/A^ in eqs. ( |60D and (|6^) (taking the real part 
and multiplying by p), respectively. 

Let us now examine in more detail the singularities in (|42|) . First note that double 
poles at half integer values do not occur. Thus the term r^^ s sin (vrp) in ( [77D is actu- 
ally not relevant. Next, to determine the residues Vp of -Bjy(z) we use the fact thatQ 
Residue {r(-2)} = (-1)^+7^!, and therefore Residue {r(-2z)} = (-l)2p+7 (2(2p)!). 
Consider the first term in (|4^ ) which is related to large-angle emission. 
p^^p^i Z/^^(— 1)^^ and there are no double poles, so 



Here r^" 



rla 



-P2p 



7rp2(2p)!r(l + /t 




In 



A2 



sin (np) + KpTT In 



A2 



Kt) 1 



COS (vrp) 



f82) 



Summing over p this gives 



large— angle 

NP 



2p 



(2p)! V Q / 



VAy^_ - ^ 1 (uA' 



n=l 



Q 



(83) 



**Since our interest is in large v values, corresponding to small t, we can replace j{—z,iy) by r(— z), 
which is a good approximation at any v > 1. 
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where we defined 
2p odd 

2p even 



^2p — — — 



A 



2p 



(3o T{l + K, 



,2 



(ingVA 

- (inQVA 



Kp sin (np) 



7T 

hi-n 1 



COS (np) 



The second term in (^), associated with colhnear emission, is treated in a similar 
manner. Note that here there are also double poles at p = 1 and 2 and consequently the 
residues have a more complicated dependence on u. The result is 



-ut\ 



col 




(85) 



where the first line corresponds to the double poles at p 



, (1) _ 

— 



Cf P2p (1 + K,p] 
Po P r(l + 



InQVA' 



the second line corresponds to the simple poles at p = 1 and 2, and the last line to the 
remaining simple poles at integers p > 3. \2p are given in 

Finally, the total non-perturbative contribution in (|67D is given by 



5Np(Q^ln v) 



,\ large- angle 
/ NP 



I NP 



^6) 



Note that in both ( [83| ) and ( pSf ) we use the same non-perturbative parameters p2p- One 
could introduce a more general parametrisation of power-corrections by defining separate 
parameters for the large-angle and collinear regions of phase-space. This, however, seems 
redundant from the current point of view: compensating the ambiguity in the perturba- 
tive sum does not require separate parameters. Moreover, the delicate relation between 
the large-angle and collinear limits of phase-space at the perturbative level suggests that 
the two are correlated in the same way also at the non-perturbative level. 

It is now straightforward to see the relation with previously suggested models for 
power-corrections 0, ^ ^ 0. As mentioned in the introduction, it was shown 
[0, ^ that the primary non-perturbative effect in the region t ^ A/Q is a shift of the 
resummed distribution to larger values of t: 



1 da 

a dt 



a dt 



(t-AiA/g). 



This is indeed the effect if one keeps only the first term in the sum in 
only non-perturbative contribution to iSnp((5^, In z/) in ( |67|) is Ai z/A/Q, and 



(87) 
): then the 



/ia = 5(t-AiA/Q). 



26 



Assuming that the average thrust is dominated by the two-jet region, the corresponding 
correction is (t) — > (t) + XiA/Q, namely it is the same parameter that controls the shift 
of the distribution and the leading power-correction to (t). This 1/Q correction to the 
average thrust was discussed extensively in the literature |21|, ^ ^ |1^, [3^ . 

\) is no longer adequate when t ~ A/Q. The 
allows parametrisation of all the terms in 



As explained in [2£, 30 1, the shift 
shape-function ^ ^ ^, 



the following identification: 

rft exp(-z/t) (Q/A) fUiQ/A) = exp 



-ut\ 
I r 



large— angle 



making 
(89) 



The single argument function {Q / A) f\^{tQ / A) is therefore a particular example of the 
two-argument non-perturbative function f^piQ'^yt) defined in eqs. (|^) and (|69D . Simi- 
larly to |30[, the moments of the large-angle shape-function am = Jo° /ia(e)e™'c^e are 
related to the parameters \2p. Writing (|89|) as 





z/eAl 




Q J 





/ia(e)cie = exp 



OO -| 

n=l 



'uAV' 



and expanding the exponents on both sides one has: 

1, cTi = Ai, o-2 = A^-A2, (Ts = A!^ - 3A1A2 + A3 



(90) 



(91) 



The function /ia(e), like the parameters A2p, depends on the regularization prescription of 
the perturbative calculation. In the absence of a non-perturbative calculation (ref. |]29[| 
gives a field theoretic definition) the shape-function and thus A„ can only be determined 
by fitting the data. 

Finally, to deal with the coUinear power-corrections, one should include (|55D. Here, 
the most important corrections arise from the terms that are enhanced by In u. Note that 
the remaining terms can be simply recast in the form 



dtexp{-ut) {QVA')a 



sub— leading (7/^2 



Win 



(92) 



similarly to the shape-function of (|89|) . These corrections become important only at 
t ~ A^/Q"^ and can be safely neglected in the region t ^ A/Q. To treat the terms that 
are enhanced by \nv we first expand the exponent. 



exp 



col 



+ 



, ^A^ /, N 

l + ^2-^(lnz/ + 7i?) + 

(In z/ ^sf - A2A4 (In p + -iEf 



A 

y (In z/ -fEf - A4 (In z/ + -/e) 



A3 
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+ 



(93) 



where we kept the sub-leading 7^ constants for convenience (see below). Next we use the 
fact that 

1 



re 

Jo 



dt- f e-"' 
t 
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to write the following convolution integral in the "+" prescription, with a generic test 
function h{t) whose inverse Laplace transform is H^u), 



Jo It 

dte-""' / rft^e^'"^* 
Jot 



h{t - 1) 
h{t - i) - h{t) 



(94) 



Now one can extract the inverse Laplace transform of the terms in (|93D, in the distribution 
sense, by expanding both sides of eq. (|9l) in powers of e. In particular, one has 



dt 
Thus, 
/coi(Q^, i) 
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(95) 



/Int 



- ln- + 0(l) 
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A^ 
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where 5'-"^ {t) is the n-th derivative of the (5-function[^ Clearly, the leading correction 
is as important around t K/Q as the large-angle corrections contained in the shape- 
function (pg]). The sub-leading terms in ( pS] ) are less important and can be ignored in a 
first approximation. 

In conclusion, the most important power correction at t 3> A/Q is the one associated 
with the simple renormalon pole at z = |. Its exponentiation amounts to a shift of 
the resummed distribution by AiA/Q. Closer to the two-jet region sub-leading power 
corrections of the form Xn{A/tQ)"' become important. Strictly based on the large /3o 
renormalon calculation, only odd values associated with renormalon poles at half integers, 
are relevant. However, since the large /Sq limit is not expected to predict the singularity 
strength of the exact Borel transform, one should expect A„(A/tQ)"' with even n to 
survive. Assuming that the exact Borel function differs from the calculated one only 
by a modification of the singularity strength, with the u dependence left unchanged, we 
predict that the even terms will be suppressed by 1/ ln((5^/A^) compared to the odd ones 
(see Allowing both odd and even terms in the sum in ( pOf ) one recovers the shape- 
function description as far as the large-angle gluon emission is concerned. Nevertheless, 
at the same time a collinear correction proportional to A2(A/tQ)^ appears in (|95|). The 
latter has a sign apposite to the relevant term in the large-angle shape-function, with 
a slightly different functional form, thus reducing significantly the overall effect of the 
non-perturbative parameter A2. 

^^^^In the last line, the appropriate integration prescription at t — > is understood. 
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Finally, we comment that an alternative way to derive the form of the power correc- 



tions from the large f3o calculation is to use the characteristic function (^41). Expanding 
the latter at small e one obtains 

^.(e) = 2v/^.+ (-i.lnQ)-3.-i.^)e+(.+ i.^ + i.^)e^/^ 

4 \ej 12 48 / V 6 72 300 / 
V36 48 480 2160 / 

+ ('_^^,5^^^6^^L^,7\ ^7/2 

V 360 3600 17640 / 

Here only non-analytic terms are associated with ambiguity of the perturbative sum (^3|), 
and therefore with power corrections One can immediately identify the terms 

{y\fe)'^ where n = 1, 3, 5, ... as the half-integer renormalon poles associated with large- 
angle soft gluons. These are the most important power corrections. As for integer 
powers e^, the only non-analytic terms are at p = 1 and p = 2, due to the logarithm. 
These ambiguities, however, are negligible for any t ^ A^/Q^ since each power of e is 
accompanied by v rather than by . Thus the picture is consistent with what we already 
learned based on the Borel analysis in the case that the singularities at integer z are just 
simple poles and are therefore cancelled by sin(7rz) factor from the analytic continuation. 

4 Relation to fixed logarithmic accuracy and fixed- 
order calculations 

In order to prove consistency of DGE with the exact exponent to NLL accuracy we will 
now compare our calculation to that of [§. This will lead to identification of coupling 



A as the "gluon bremsstrahlung" coupling |T^. Next we briefly recall the way to match 



the resummed result to the fixed-order calculation. 



4.1 NLL resummation from an evolution equation 

The calculation of the differential cross-section by DGE requires to identify the cou- 
pling A. The analogy with the Abelian theory |l^, |14|, |16| suggests that there is a unique 
all-order definition of A (analogous to the Gell-Mann-Low effective charge) which does 
not depend on the observable considered, such that @ is the leading term in the skeleton 
expansion. Since a systematic skeleton expansion has not yet been constructed in QCD 
other considerations must be used. The analogy with the Abelian theory allows one to 
fix the large /5o coupling, as done in (^). To go beyond this limit we require consistency of 
the exponent to NLL accuracy. We will see that this amounts to identifying the coupling 
A as the "gluon bremsstrahlung" coupling []T^. 

The standard way of calculating the logarithmically enhanced cross-section is based 
on an evolution equation for the jet mass distribution J((5^, A;^), where the kernel is the 
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splitting function of a gluon off a quark, P(a, x). At leading order the splitting function 
is 



P{a,x) = Cpd 



Cf a 



X 



2(1+^) 



(96) 



where a{k'^) is the coupling {as{k'^)/Ti) and x is the longitudinal momentum fraction in 
the branching. Fixing the coupling to be the "gluon bremsstrahlung" coupling ||17|| , 



:/5o + 



3 ~ 12 



(97) 



this splitting function is correct to next-to- leading order, as far as the singular 1/(1 — x) 
terms are concerned. At NLL accuracy, the evolution equation is H], 



c/lnQ2 



dxP (S((l - x)Q^),x) [e-''(^-") - 1] J,(g2) 



(9J 



where JuiQ"^) is the Laplace transform of J{Q^, k"^) with respect to k"^ jQ"^. The solution 
of this differential equation is 



lnJ,(g' 
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(99) 



To this accuracy 1 — thrust, t, is the sum of invariant masses of the two hemispheres, so 



/■CO _ _ , _ 

i?(Q2^ t) = dedk^J{Q\ k^)J{Q\ ¥)d (tQ^ 



and therefore, 



du 



■e^'^exp {2 In J^(g2) I 



(100) 



(101) 



res Jc 27^1/ 

From here, a straightforward calculation gives the leading and next-to-leading logs in the 
thrust distribution. 

Let us now compare this calculation with ours. Specifically, 2ln Ji^lQ"^) in (|101|) 
should be compared with in (^^. The two quantities are defined in (^) and 

(|35D, respectively. The first observation is that the integration variable u in (|99|) can be 
identified with t, making the two similar provided one identifies 



1 da 
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or, explicitly (see |29|) , 



*-^eff(6g^) 

t2 e 



2 _ 1 _ f! 
t t2 ^3 



t Jo 



2 



t Jo 
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X 



(102) 



1 — X 



2(1 + ^) 



(103) 

If one identifies y4eff(eg^) with /3oa((l — x)tQ'^), and changes variables e = (1 — x)t, the 
r.h.s. becomes 



Jt'2 e 



2 _ 2e 
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One thus finds that the first term in the square brackets on both sides is the same. Recall 
that this term is responsible for the leading logs. The other terms, contributing to sub- 
leading logs, are not the same. However, the NLL are: to NLL accuracy it is enough to 
consider the conformal part |^ by fixing the scale of the coupling. The latter gives, in 
both cases, — j|y4((5^) from the collinear limit e = t. 

We thus learned that DGE is consistent with the exact exponent to NLL accuracy pro- 



vided one identifies the running coupling with the "gluon bremsstrahlung" coupling |jT^ 
of eq. (P?!). In evaluating the Borel sum ( |55D we therefore use 

Mk') = . r. . ,.1 (104) 



instead of (g). This simply amounts to fixing the definition of A (see 

The qualitative features of the SDG characteristic-function based calculation appear 
also in the splitting-function based calculation: there appear e/t^ and e^/t^ terms which 
contribute to the log-enhanced cross-section due to the collinear limit e ~ t. As we 
saw in the previous sections these terms lead to factorial enhancement of sub-leading 
logs (already prior to exponentiation). However, the actual values of sub-leading logs 
predicted by the splitting-function beyond the NLL are different. Performing the integral 
on the r.h.s. of (|103|) explicitly we arrive at the following expansion for ^^(Q^,^) prior 
to exponentiation: 

-^{Q\ t) = % 1(2. L - 1.5) A{Q^) + {3L^ + 1.8333 L - 4.2500) A{Q^f (105) 
a dt t Po 

+ (4.6667 + 8.5000 - 9.5242 L - 8.8151) A{Q'^f 

+ (7.5000 + 21.833 - 17.359 - 50.085 L - 9.957) A{Q'^f 

+ (12.400 + 48.500 - 31.339 - 194.72 - 95.096 L + 10.97) A{Q'^f + ■■■ 

This expansion can be compared directly with (|25|) . The LL and NLL are of course 
the same. Sub-leading logs are similar in many cases; note in particular the NNLL. 
Nevertheless, at lower logs there are big variations, as one can deduced directly from ([1031) . 



The fact that the two calculations differ beyond NLL is of no surprise. Most im- 
portantly, the kinematic approximations made in the splitting-function based calculation 
(e.g. the matrix element, the definition of the observable), modify the logs beyond NLL 
accuracy. 

It should be emphasised that some approximations have been made also in the 
characteristic-function based calculation. Although this calculation is exact at the level 
of three partons with an off-shell gluon, there are approximations beyond this level: the 
massive-gluon characteristic-function does not take into account the possible branching 
of the gluon into opposite hemispheres, the so-called non-inclusive contribution. When 
the gluon is roughly collinear with one of the quarks (the only limit contributing to sub- 
leading logs in ^^{Q'^,'t) ) this type of branching is probably not significant. Indeed, 

the comparison between ^^{Q'^,t) and the leading term in Pq in the (numerically 



SDG 
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evaluated) exact NLO coefficient in sec. 2.2 con&ms this assertion. Another important 
aspect is the vahdity of the exponentiation formula (|37| ) beyond NLL accuracy. As dis- 
cussed in sec. 3, correlations in the emission of multiple soft and coUinear gluons may 
influence the sub-leading logs. This effect was neglected. 



NLO 
res 



4.2 Matching to the fixed-order calculation 

At large values oft = 1 — thrust, the resummed cross-section (be it DGE or NLL) must be 
complemented by the exact fixed order result. To do this we will use the log-R matching 
scheme. The basic idea of matching is to replace the low orders in the resummed expres- 
sion by the exact low order terms. Schematically one writes -Rmatched = R^^° + Rres — Rj 
where i?res is the resummed cross-section and R^es° expansion to NLO. There is no 
general argument for matching the integrated cross-section itself (R-matching) as com- 
pared to matching some function of R. Exponentiation may hold beyond the level of 
the logarithmically enhanced cross-section, so it is natural to match the logarithm of R 
(log-R matching). In practice the difference between different matching procedures is 
small (see e.g. 0) and for simplicity we only consider log-R matching. 
Using the known NLO calculation [^, [1^ , 

R{Q^t) = 1 + Ri{t)a + R2{t)a^ + ■ ■ ■ 



where Ri{t) is known analytically and R2{t) is parametrised based on a numerical com- 
putation by EVENT2 Il3|, we write 



R{Q\ t) = exp ^RMQ ) + [R2 - ^Ri) a(Q ) + In RiQ\ t) - In R{Q', t) 

(106) 

Here In RlQ"^ ,t)\^^ is our resummed result for the logarithmically enhanced integrated 
cross-section, given by eqs. (|6^) and (|63| ) and \nR{Q'^,t)\^^^ is its expansion, up to NLO 
in MS, 
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The matching (|106| ) is done in the MS scheme. 



5 The phenomenological implication 

We saw that a careful treatment of running coupling effects (or renormalons) in the 
logarithmically enhanced cross-section requires resummation of sub-leading logs. Our 
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investigation in section 3.2 leaves no doubt: sub-leading logs contribute significantly 
at all physically interesting values of the thrust. On the other hand, it is well known 



that successful fits to the distribution were obtained in & 30, 29 with a much sim 



pier treatment of the perturbative calculation, based on the standard NLL resummation. 
Non-perturbative corrections whose magnitude is not under control theoretically are any- 
way important, and must be included in any fit. One may argue that a more complete 
perturbative treatment is not necessary, and that it amounts simply to a redefinition of 
the non-perturbative correction. The answer is, of course, that as soon as quantitative 
predictions are expected one should start with a reliable perturbative calculation. Phys- 
ically, the transition between the perturbative and non-perturbative regimes is smooth, 
and one must make sure that an arbitrary separation of the two does not have too strong 
an impact on the results. Legitimate separations are cut-off regularization of momentum 
integrals and principal- value Borel summation [Q. It was shown in that these two 
procedures can be fully consistent with each other provided that renormalon resummation 
is performed. This means in particular that the extracted value of as is the same using 
the two prescriptions. Such consistency cannot be achieved if a truncated perturbative 
expansion, e.g. NLO, is used instead. Our results in sec. 3.2 imply that NLL truncation 
of the exponent is just as dangerous. By performing renormalon resummation, such a 
truncation is avoided and a more reliable comparison with data can be achieved. 

Here we shall perform global fits to the thrust distribution at all energies, based on our 
perturbative calculation of section 3. As discussed in section 3.3 we use Borel summation. 
An equivalent calculation can be done using cut-off regularization. The translation of 
our results to this language is straightforward and may be useful once comparison can 
be made to a non-perturbative calculation of the shape-function with a hard cut-off as 
an ultraviolet regulator. 

The perturbative calculation must be supplemented by non-perturbative corrections 
in order to be compared with the data. Our guiding principle is that the functional 
form of the non-perturbative corrections is dominated by the ambiguities identified in 
the calculation of the perturbative sum. We assume that the SDG based calculation is 
sufficient to identify these ambiguities. Of course, there are inherent limitations to this 
calculation: splitting of a gluon into partons that end up in opposite hemispheres is not 
taken into account, nor are correlations between the emitted gluons. These effects are 
not expected to be important in the two-jet region, so it is reasonable to neglect them. 



5.1 Perturbative calculation 

Before attempting any fits to the data involving non-perturbative corrections, it is use- 
ful to examine the impact of the resummation simply by looking at the perturbatively 
calculated cross-section in fig. |[ Needless to say, the fixed-order calculations, LO and 
NLO, do not describe the physical cross section in the two-jet region. The large differ- 
ence between the NLO and the LO signals large higher-order corrections. So does the 
renormalization scale dependence of the result (not shown in the plot). The qualitatively 
different behaviour of the NLL result is clearly seen: thanks to exponentiation it van- 
ishes at t = 0. Nevertheless, as fig. ^ shows, the renormalization scale dependence is 
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Figure 4: DGE, the standard NLL and the LO and NLO results as a function of t in the 
two-jet region, at Q = Mz. Fixed-order and NLL results are in the MS scheme. As an 
example of scale dependence the NLL result is shown at two different renormalization 
points fJ^R = Q and fiR = Q/2. For the LO and NLO results = Q. We assume 
af^ (Mz) = 0.110. K data is shown for orientation. 

still appreciable^]. We stress that the factor 2 between the renormalization points is com- 
pletely arbitrary. This difference should not be considered as a reliable estimate of the 
error due to running coupling effects that have been neglected, since the physical scale, 
(e.g. the transverse momentum) in the two-jet region is significantly smaller than Q. As 
anticipated in sec. 3.2, the impact of the additional resummation of sub-leading logs by 
DGE is quite significant. Having performed this resummation, the renormalization scale 
dependence (at the level of the logarithmically enhanced cross-section) is avoided. 

In fig. |, as in all forthcoming calculations we use the so-called log-R matching 
scheme , defined in Eq. ( |106D , to complement the calculation of the logarithmically 



*It has been shown |T6j that renormalon resummation in single scale observables can be imitated 
by an appropriate choice of the renormalization scale, similarly to the original formulation of the BLM 
approach [|l5j. One might wonder to what extent DGE can be imitated by simply tuning the scale of 
the coupling in the standard NLL formula. Physically, the naive scale-setting approach is bound to fail 
for multi-scale observables such as differential cross-sections, since the typical scale of emitted gluons is 
not anymore a fixed fraction of the center-of-mass energy, but rather some function of all the external 
parameters. In our case, it is a function of both Q and t. Indeed, choosing lower scales for the coupling 
in the NLL formula one cannot reproduce the DGE result. 
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enhanced cross-section. Renormalization scale dependence appears only beyond the level 
of the logarithmically enhanced cross-section and it is therefore suppressed by t. Numeri- 
cally, the residual scale dependence (see table [|) is very small as long as the log-enhanced 
cross-section dominates. The difference is much too small to be observed in the resolution 
of fig. ^. At large t {t ^ 0.2) the terms which are suppressed by t start making an impact 
and eventually the scale dependence becomes comparable to its values in the NLL and 
the NLO results. 



Table 1: Residual renormalization scale dependence: the relative difference between the 
differential cross-section at Q = Mz as calculated from ( |106D when the matching with 
the NLO result is performed at fiji = Q and at Q/2. 



t 0.02 0.05 0.10 0.15 0.20 0.25 0.30 0.33 
% 0.20 0.15 0.01 0.29 0.58 1.1 2.8 7.3 



Next, we confront the question of reliability of the approximation technique we used in 
section 3.3 in evaluating the Borel sum. In the calculation presented in fig. ^ we evaluated 
the Borel sum in (^) using Pmax = 8- This means that the function V{z) in eq. (^) was 
approximated by 8 poles. Note that pmax is also the largest power of In z/ respected by 
this approximation (see (^)). In addition, the renormalon integrals /^(InQ^/A^) were 
evaluated using the two-loop running coupling in the "gluon bremsstrahlung" scheme 
(see section 4). How sensitive are the results to these approximations? To answer this 
question we repeat the calculation varying the choice of pmax and the running coupling. 
Considering first the choice of Pmax, we show in fig. ^ approximants of increasing order. 
The convergence is quite remarkable. In the upper plot the Pmax = 5 curve (dashed) can 
hardly be distinguished from the Pmax = 8 curve. Indeed, as the lower plot shows, the 
difference between them is very small and it exceeds 1% only around the distribution 
peak. For our fits, we choose Pmax = 8 which guarantees a truncation erroiQ of less than 
1% down to t = 0.01 at Mz (or 0.01 Mz/Q at Q). 

The stability of the calculation with respect to the running coupling can be checked 
by repeating the calculation with a different renormalization-group equation. In prin- 
ciple, one may also consider changing the definition of A, thus replacing the "gluon 
bremsstrahlung" coupling by some other coupling that matches the same Abelian limit. 
This was done in a similar context in the case of the average thrust . It was found there 
that such a replacement introduces a non-negligible variation in the calculated observable, 
but the final impact on the extracted value of a.^ was rather small (of order 1%). In the 
context of the distribution, the choice of the "gluon bremsstrahlung" coupling is almost 
imposed on us by consistency with the NLL exponentiation kernel. Therefore we shall 
not make such a modification. The simplest variation in the renormalization-group equa- 
tion is to replace it by the one-loop. In this case the renormalon integrals /p(lnQ^/A^) 
are evaluated by eq. (pOD instead of eq. (^). Such a comparison is presented in fig. ||. 
The difference between the two curves is quite significant, in particular in the peak re- 

^This truncation error must not be confused with the ambiguity of the perturbative sum. 
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Figure 5: The convergence of increasing approximants to the perturbative cross-section, 
according to the technique of section 3.3, eq. ([631). The upper plot shows the differential 
cross-section for Pmax = 3,5 and 8 and the lower plot shows the relative error of the 
differential cross-section with respect to the Pmax = 8 calculation, for pmax = 5,6 and 7. 
We fix af^ (Mz) = 0.110 and use two-loop running coupling; log-R matching is applied. 



gion. Note, however, that contrary to the two-loop case, the one-loop resummed result 
does not include all the next-to-leading logs (the terms proportional to Pi are omitted). 
The difference between two- and three-loopf] running coupling is expected to be much 
smaller. However, it may well be non-negligible close to the peak, so it is certainly worth 
consideration in the future. In the following we will use the two-loop running coupling. 



5.2 Non-perturbative corrections 

It was realized before [0, |2^, |2B|, ^ that non-perturbative power corrections of the form 
An(A/tQ)" to any power n become relevant around the distribution peak. As shown in 
section 3.5, such corrections are indeed expected based on the form of the ambiguity of 
the perturbative result. If all the powers are as important, one would need an infinite 
number of parameters - or a completely arbitrary shape-function - to bridge the gap 
between the perturbative result and the data. In practice one assumes that there is a 



Hn the spirit of and our discussion in sec. 4, the appropriate value for P2 should probably be set 
based on the NNLO coefficient of the singular term 1/(1 — x) in the splitting function. 
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Figure 6: The resummed cross-section at Mz based on one- loop (dashed) and two-loop 
(full line) running coupling. We fix af^ (Mz) = 0.110 and use pmax = 8; log-R matching 
is applied. K data is shown for orientation. 



region in t where a power corrections expansion such as (0) is dominated by its first 
few terms, or in other words, that the shape-function can be described by its first few 
moments. It is therefore essential that An do not grow too fast with n. Having made this 
assumption, one can put a cut on the fitted range in t for any given energy, and fit the 
data for t > ti using a reasonably fiexible parametrisation of the shape-function. This 
cut should scale like the distribution peak, i.e. as l/Q, 

tL{Q) = tL{Mz) (108) 

In any case, the dependence of the results on the precise location of this cut and on the 
specific functional form of the shape-function must not be too large. 

The simplest possibility is to place the cut at large enough values of t such that the 
only relevant non-perturbative parameter will be Ai, which controls the shift of the distri- 
bution (pTl). This approach was suggested in 0. The advantage is that the perturbative 
calculation is better under control and the "model dependence" in the parametrisation of 
non-perturbative corrections is minimized. The price, however, is that one does not use 
a large part of the data, and in particular, one excludes the region where the distribution 
changes sharply and is therefore expected to be most constraining and most informative. 
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Here we use various different minimal t cuts and perform fits where a flexible shape- 
function is convoluted with the perturbative distribution as well as fits with just a shift. 
For the shape-function-based fits the minimal cut can be quite low, thus using almost all 
available data. At the low end we choose ^^.(Mz) = 0.01. This is equivalent to Qt ^ 4A 
and it guarantees convergence of the Borel function approximation within 1%, as shown 
in fig. ^. Less ambitious cuts ^^.(Mz) = 0.02 and 0.03 are useful to check the stability, 
which is indicative of whether the power-correction expansion is convergent enough. For 
the shift-based fit higher cuts must be used. Also here stability of the results is important 
to guarantee that sub-leading power corrections are not too relevant in the fitted range. 

For simplicity we shall assume that the parameters of the shift or the shape-function 
are Q independent, even though we have shown that logarithmic dependence is expected 
in general. In principle, a fit that allows logarithmic dependence of A„ according to 
eq. ( p4D is more appropriate. But, since Kp are not known, the price will be having 
more parameters which can only be constrained by very precise data at several different 
energies. The available data is not constraining enough for such a fit. 

5.2.1 Fits with a shifted distribution 

Let us consider first a fit based on shifting the perturbative distribution. We use ex- 
perimental data from all energies in the range 14 GeV < Q < 189 GeV (summarised in 
table ^ below). When available, the systematic errors have been added in quadrature to 
the statistical errors. Table § summarizes the fit results for and Ai as a function of 
the low cut. 

Table 2: Fits of as and Ai to data in the range 14 GeV < Q < 189 GeV with a lower 
limit tiiMz) Mz/Q and upper limit tn = 0.32. The last two columns show the obtained 
value of P^r degree of freedom and the total number of data points used in the fit, 
respectively. 



tL(Mz) 






^(Mz) 


Ai 


xVdof 


points 


0.02 


0, 


,1088 


± 0.0002 


2.71 ± 0.11 


0.97 


231 


0.03 


0, 


,1081 


± 0.0003 


3.11 ± 0.16 


0.80 


213 


0.04 


0, 


,1085 


± 0.0004 


2.81 ± 0.23 


0.68 


193 


0.05 


0, 


,1092 


± 0.0005 


2.32 ± 0.30 


0.66 


177 


0.06 


0, 


,1100 


± 0.0006 


1.76 ± 0.40 


0.64 


165 


0.07 


0, 


,1106 


± 0.0008 


1.39 ± 0.48 


0.65 


152 


0.08 


0. 


,1108 


± 0.0009 


1.27 ± 0.53 


0.62 


140 


0.09 


0. 


,1114 


± 0.0012 


0.84 ± 0.69 


0.62 


127 


0.10 


0, 


,1119 


± 0.0014 


0.57 ± 0.81 


0.64 


119 


0.11 


0, 


,1126 


± 0.0019 


0.18 ± 1.09 


0.64 


109 


0.12 


0, 


,1130 


± 0.0029 


-0.08 ± 1.26 


0.61 


101 


0.13 


0. 


,1156 


± 0.0037 


-1.26 ± 1.67 


0.60 


94 


0.14 


0. 


,1123 


± 0.0045 


0.32 ± 2.48 


0.60 


91 
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Figure 7: The correlation between Og and Ai, as determined by a shift-based fit with 
various cuts (full one-sigma contours) and by a shape-function-based fit with a cut 
^^(Mz) = 0.01 (dashed one-sigma contour). The additional lines indicate the values of Ai 
in fits with a fixed Og- The full line corresponds to the shift-based fit with t^lMz) = 0.05 
with and the dashed to the shape-function-based fit with ^^(Mz) = 0.01. When is 
free (the centre of the corresponding contour) x^/dof = 0.66 and 0.84, respectively. The 
points where x^/dof = 1 are indicated. 

First of all we note that the fits are very good: the x^/dof values are unusually low. 
Next, we see that the results depend on the low cut. Particularly stable results, as far 
as the extracted value of is concerned, are obtained for t^^Mz) ^ 0.06. Above this 
value consecutive cuts yield ag values that are within errors of each other. This stability 
means that sub-leading power corrections in this region are small enough to be neglected. 
Thus we regard the result af^ (Mz) — 0.110 as rehable. Too high cuts tL(Mz) > 0.10 
are disfavoured because the number of fitted points decreases significantly (note that the 
propagated experimental error increases accordingly). Cuts hke ^^(Mz) = 0.05 have the 
advantage of increased sensitivity to Ai. Lower cuts are probably too sensitive to sub- 
leading power corrections. This is reflected in the stronger dependence on the cut and in 
the deterioration of x^/dof. In the following we shall use ^^(Mz) = 0.05 for illustrations. 
It has the additional advantage of consistency with the shape-function-based fit (see 
below) . 

The table and fig. |^ show that there is a strong correlation between ag and Ai. A 
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Figure 8: Comparison with experimental data at Q = 14, 22, 35 and 44 GeV with the 
perturbative calculation (dotted) and the corresponding fits using a shift (dashed) or a 
shape-function (full line). The perturbative calculation (DGE) is performed using Borel 
summation with two-loop runing coupling. The fits are based on the range ^^(Mz) ^ < 
t < tff = 0.32, where t^iMz) = 0.05 in the case of a shift and ^^(Mz) = 0.01 in the case 
of the shape-function. The coupling is af^ (Mz) = 0.1093 (0.1092) which is the best fit 
value in the shape-function (shift) based fit. 

similar correlation was found in 0. This means that quantitative discussion of the values 
of non-perturbative parameters, like is relevant only for a fixed a^. 

In all the fits discussed so far we used an upper cut tn = 0.32 which is Q independent. 
We checked that the location of this cut does not infiuence the fit. For example, keeping 
the lower cut fixed ^^.(Mz) = 0.05 and changing the upper cut between tn = 0.29 and 
tff = 0.35 (corresponding to a total of 154 to 184 fitted points, respectively) the central 
value of as changes by just 0.2% (and that of Ai by 4%). 

Figures ^ through |T0| show the fitted distribution and the original perturbative distri- 
bution together with the experimental data. For large Q we show only the small t region, 
which is more interesting. The agreement between the shifted distribution and the data 
is good. Nevertheless, it is clear that around the peak additional corrections are required 
(in the case of Mz, for example). Indeed we saw that such corrections are expected based 
on the renormalon ambiguity. 
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Figure 9: Comparison with experimental data at Q = 91 and 133 GeV with the per- 
turbative calculation (dotted) and the corresponding fits using a shift (dashed) or a 
shape-function (full line). See fig. || for further details. 

5.2.2 Fits with a shape-function 

Consider now a fit where the perturbative distribution is convoluted with a shape- 
function. According to sec. 3.5, this means, in fact, two convolutions: one with the 
large-angle shape-function /ia(e), defined in eqs. ( p9|) and (|90|) , and another which cor- 
rects for the collinear contribution where we only include the leading term in 



fUQ^t) = Sit)-\2 



QH 



5«(t). 



The main difference with previous works where a shape-function was used [^y], ^ is in 
the starting point: the perturbative distribution there is the standard NLL result (with 
some cut-off regularization) while we start with the DGE result avoiding the truncation 
of sub-leading logs. The latter provides an analytic continuation of the perturbative 
treatment beyond its original range of applicability. 

In spite of the significant difference in the perturbative distribution as compared 
to |30|, we consider their ansatz for the shape-function appropriate also in our case. 



barring the fact that it is positive definite. To allow in principle the suppression of even 
moments as implied by our analysis in sec. 3.5, we consider a function that is flexible 
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Figure 10: Comparison with experimental data at Q = 161, 172, 183 and 189 GeV with 
the perturbative calculation (dotted) and the corresponding fits using a shift (dashed) or 
a shape- function (full line). See fig. ||for further details. 



enough: 

/ia(e) = noe'^ (1 + he + he^) e-^^^^^^' (109) 

where no is fixed by the normalization requirement / f{t)dt = 1 and q,ki,k2,bi and 62 
are free parameters to be fixed by the fit to the data. We assume g > so that /ia(0) = 
and 62 > (or 62 = and bi > 0) so that the shape-function decreases at large e. 

Note that with a shape-function of this form, the central moments A„ do not have a 
definite sign. In particular both positive and negative values for Ai are possible. When 
a cut-off regularization is used it is intuitively clear that the shift is to the right: soft 
gluons widen the jets. Indeed Ai was found to be positive in both the shift-based fit in 
P] and in the shape-function PSI, RUl based fit. The shift parameter Ai was shown to 



be roughly consistent with the leading power correction to average event-shapes p, p9| , 
which is always positive. The intuitive argument is not anymore valid when principal 



value regularization is used instead. According to in this regularization Ai is still 
positive. We stress, however, that there is no general reason to exclude a priori the 
possibility of a negative Ai. 

Performing a fit with /ia(e) and the collinear correction ( |109| ) with a lower cut t/^(Mz) = 
0.01 and an upper cut iff = 0.32, we obtained a very good fit with x^/dof = 0.84 (for a 
total of 246 fitted data points). The fit results are shown in figures | through |1^. The 
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Figure 11: Comparison of DGE convoluted with the best fit shape-function (62 = 0) 
with experimental data for Q = 14 up to 189 GeV. Both the DGE result and the data 
have been multiplied with a factor 20^, with = 0..10 for Q = 14. .189 GeV. Here 
af^ (Mz) = 0.1093. Data points with large errors (> 100%) are not shown. 

entire range of t is shown for all energies in fig. The parameters in this case are 

af{Mz) = 0.1093 ± 
q = 0.580 ±0.115 
bi = 0.408 ±0.045 
&2 = 

corresponding to 

Ai = 2.34 ±0.21 (111) 
A2 = 0.37 ±1.38. 

It turns out that the fit has a strong preference for a linear term in the exponent, so that if 
both bi and 62 are free, then 62 is set by the fit to zero. In order to try another functional 
form, thus testing the sensitivity of the fit results to the properties of the shape-function, 
we also fitted the data in the same range with the function where the linear term in the 
exponent is missing: 61 = 0. In this case we got a best fit with x^/dof = 0.88, i.e. shghtly 
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0.0004 



= -0.099 ± 0.034 
A;2 = 2.5 ■ 10-6 ± 0.003 



(110) 



worse than the previous one, with the following parameters, 



af (Mz) = 0.1082 ±0.0003 
q = 1.53 ±1.23 

61 = 

62 = 0.099 ± 0.058 



ki = -0.422 ± 0.052 
k2 = 0.052 ±0.010 



(112) 



corresponding to 



Ai 
A2 



2.92 ±0.14 



(113) 



-3.52 ±0.57. 



We see that most of the parameters of the shape-function are fairly well constrained 
by the fit. Unfortunately, this does not imply that the central moments of the shape- 
function are well determined. As explained above, due to the correlation between as and 
the non-perturbative parameters, it makes sense to fix «<, in the fit when comparing the 
results for the two possible exponents. Doing so with = 0.1093 we still obtain good fits 
in both cases. For a linear term in the exponent we get Ai = 2.35 ± 0.05 (x^/dof = 0.84) 
and for a square in the exponent we get Ai = 2.45 ± 0.04 (x^/dof = 0.92). Now the 
values of Ai are much closer. Note that also the propagated errors become smaller, since 
as and Ai are so strongly correlated. 

Figure shows how Ai varies when the fits with the shift and with the shape-function 
(with a linear term in the exponent) are performed with a fixed value of as- As intuitively 
expected, forcing a larger value of as, Xi decreases. For af^ (Mz) — 0.109 similar values 
are obtained for Ai from the shift and from the shape-function (note that this depends 
on the chosen lower cut), however, for higher values of the coupling the values of Ai start 
differing. In this case the shape-function is not so well approximated by a (5-function (a 
shift) and its higher central moments A„ play a role. 

To allow a quantitative discussion concerning the higher central moments let us fix 
also Ai and compare the fits with the two exponents for the same as and Ai. We choose 
them to be the values obtained in our best bit (|TToD : af^(Mz) = 0.1093 and Ai = 2.34. 

Quantitative information on higher moments is theoretically important for several 
reasons. First, it would be interesting to see the effect of the collinear contribution. 
At the same time it would be interesting to check our prediction that A2 (like other 
even central moments of the shape-function) is suppressed. Even ignoring the strong 
dependence on as, these two ingredients are hard to address by fitting the data: if A2 
is small then the collinear correction is not important, and conversely - if the collinear 
correction is included a major part of the effect of A2 is cancelled, so A2 would be very 
hard to access. To deal with this we make three different fits with the same fixed values 
of as and Ai mentioned above: with the collinear correction, without it, and finally a fit 
where A2 is fixed to zero. The results are summarised in table 

The first conclusion is that A2 is not well determined by the fits: even for fixed values 
of as and Ai the two functional forms yield somewhat different results. Nevertheless, a 
vanishing A2 is certainly not inconsistent with the data: the fit does not deteriorate much 
by fixing it to zero. One also sees that the collinear correction is not important: including 
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Table 3: Determination of A2 and A3 for fixed af^ (Mz) = 0.1093 and Ai = 2.34 





62 = 


h 


= 








A2 


A3 


xVdof 


A2 


A3 




xVdof 


with collinear 


0.42 ±0.54 


-52± 12 


0. 


84 


-0.75 ±0.43 


-20 ± 


7 


0.95 


no collinear 


0.48 ±0.52 


-53± 12 


0. 


84 


-0.69 ±0.42 


-21 ± 


7 


0.95 


A2 = 





-44 ± 5 


0. 


84 





-31 ± 


3 


0.96 



it or not leaves A2 almost identically the same. Had the value of A2 been significantly 
different from zero, such a situation would have been impossible, but for a small (and 
not very well constrained) A2 it is quite reasonable. 

Table ^ also gives the A3 values. It is clear that A3 is not constrained enough unless 
as, Ai and A2 are all fixed. Still, one thing that can be said in general (varying the other 
parameters, including as) is that negative values for A3 are strongly preferred. Thus, 
a positive-definite ansatz for the shape-function (in the principal-value regular isat ion!) 
is not adequate. Note that the preferred signs of Ai (positive) and A3 (negative) are 
compatible with p2p being all positive (see eq. (|^) and (Q)). 

The shape-function-based fit can be considered an alternative way to extract the 
value of as- We saw already, that in spite of the fact that the determination of the 
non-perturbative parameters is far from being complete, as does not change much. In 
particular, for the two functions used, af^ (Mz) varies only from 0.1082 to 0.1093. A 
further stability check is to what extent it is sensitive to the location of the low cut 
tiiQ)- The answer is provided by table ^. The sensitivity to in the determination 

Table 4: Fits of as and the shape-function (with a linear term in the exponent, 62 = 0) 
as a function of the lower cut tiiQ). The upper cut is tn = 0.32. 



tL(Mz) 


«r (Mz) 


Ai 


A2 


xVdof 


points 


0.010 


0.1093 ± 0.0004 


2.34 ± 0.21 


0.37 ± 1.38 


0.84 


246 


0.020 


0.1101 ± 0.0005 


2.06 ± 0.22 


2.61 ± 1.56 


0.73 


231 


0.030 


0.1109 ± 0.0007 


1.84 ± 0.25 


3.41 ± 1.71 


0.70 


213 


0.040 


0.1108 ± 0.0011 


1.88 ± 0.31 


2.80 ± 1.84 


0.63 


193 


0.050 


0.1102 ± 0.0012 


2.04 ± 0.45 


4.23 ± 2.56 


0.63 


177 



of as is rather small. In particular it is smaller than in the case of a shift-based fit. 
There is still some tendency that as becomes larger for higher (thus less ambitious) 
cuts. While higher cuts are safer from the point of view suppressing sub-leading power 
corrections, they are less constraining as far as the multi-parameter shape-function is 
concerned. Note also that the propagated experimental error on as and Ai increases 
for higher cuts. "Ambitious cuts" where the peak region itself is included in the fit , 
e.g. ^^(Mz) = 0.010, are therefore favoured. We conclude that according to the shape- 
function-based fit af^ (Mz) ranges between 0.108 and 0.111. These variations should be 
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considered part of the theoretical uncertainty. 

Finally, we summarise in table § the contribution to in the fit by each experiment 
separately. The middle columns correspond to a fit of the shifted distribution for as and 
Ai with a lower cut ^^(Mz) = 0.05 and an upper cut tu = 0.32. The columns on the 
right correspond to a fit for Ug and the parameters of the shape-function with a lower 
cut tL(Mz) = 0.01 and upper limit tn = 0.32. The OPAL and SLD data sets at 91.2 
GeV are not included in the fit since they are not consistent with the other data sets at 
this energy for small t. We note that there are no particular trends in the x^/point as a 
function of energy. 



Table 5: Contribution to the two fitting procedures. 



Experiment Reference Q [GeV] 


(shift) 


points(shift) 


X' (SF) 


points (SF) 


TASSO 


53 


14.0 


0.00 





4.58 


7 


TASSO 


53 


22.0 


1.43 


2 


3.17 


8 


JADE 


54 


35.0 


3.80 


6 


29.31 


11 


TASSO 


51 35.0 


1.54 


4 


4.10 


9 


JADE 


54 


44.0 


5.57 


7 


14.47 


11 


TASSO 


53 


44.0 


9.11 


5 


12.86 


9 


AMY 


55 


55.0 


2.49 


4 


10.75 


5 


ALEPH 


56| 91.2 


4.03 


10 


8.12 


17 


DELPHI 


57 


91.2 


9.68 


13 


17.18 


17 


L3 


58 


91.2 


5.74 


7 


6.13 


9 


MARKII 


59 


91.2 


7.29 


6 


9.54 


9 


OPAL 


60 


91.2 


(27.07) 


(27) 


(286.09) 


(31) 


SLD 


61 


91.2 


(9.10) 


(6) 


(96.14) 


(8) 


ALEPH 


62 


133.0 


6.08 


7 


6.28 


8 


DELPHI 


63| 133.0 


5.80 


7 


6.02 


8 


OPAL 


64 


133.0 


2.66 


7 


3.44 


10 


ALEPH 


65 


161.0 


5.62 


7 


12.34 


8 


DELPHI 


63| 161.0 


9.96 


7 


12.19 


8 


OPAL 


66 


161.0 


2.64 


8 


3.41 


10 


ALEPH 


65 


172.0 


5.94 


7 


6.22 


8 


DELPHI 


63 


172.0 


5.08 


7 


5.58 


8 


L3 


67 


172.0 


3.81 


10 


4.14 


11 


OPAL 


68 


172.0 


3.50 


8 


3.77 


10 


ALEPH 


69| 183.0 


4.53 


7 


8.35 


8 


DELPHI 


63| 183.0 


7.54 


15 


8.40 


17 


OPAL 


68 


183.0 


0.97 


8 


1.45 


10 


OPAL 


68 


189.0 


0.96 


8 


1.68 


10 


Sum 


115.8 


177 


203.5 


246 
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5.3 Moments of the thrust distribution 



Given the assumption that the first few moments of the thrust (t") are dominated by the 
two-jet configuration, and therefore by the logarithmically enhanced cross-section, one 
can relate ^ the parameters to the power corrections of the moments: 

(^)two-jet ^ (t)pT + AlA/g (114) 

I two— JCt \ / PT - ' ±-1 ' 



These predictions have been compared in |32] with the power correction from a single 



gluon emission, based on the characteristic function (|T^, 

(i)soG ^ (t)p, + AA/Q + 0(AVQ=^) (115) 
t^) - if) +0{X'/Q''). 

I SDG \ / PT ^ / -V / 

The assumption that the two-jet region is dominant probably holds for the average thrust. 
The immediate conclusion is that the shift of the distribution Ai must coincide with the 



leading power-correction for the average, A. The latter was determined from a fit in 114 
using the principal value Borel sum regularization. For a^^ (Mz) = 0.110, which is the 
best fit value in the case of the average thrust, AA = 0.62 ± 0.12 GeV. This number can 
be compared with our current fits to the distribution with a fixed af^ (Mz) = 0.110. We 
get AiA = 0.44 ± 0.02 GeV in a shift-based fit with a cut tL(Mz) = 0.05 (corresponding 
to x^/dof = 0.67 for 177 points) and AiA = 0.48 ±0.01 GeV with a shape-function-based 
fit with a cut ^^(Mz) = 0.01 (corresponding to x^/dof = 0.85 for 246 points). The power 
correction extracted from the distribution is somewhat lower than the corresponding 
power correction extracted from (t), however the discrepancy is not large. Note that 
there is a difference between the resummation procedures used in the two cases. In the 
case of the distribution DGE is used, where the exponentiation plays a major role but 
terms that are suppressed by t are neglected. On the other hand in case of the average 
a single dressed gluon renormalon sum is used, where terms that are suppressed by t 
(and higher powers) are included but multiple gluon emission is neglected (apart from 
the NLO term). 

An even more difficult case is that of the second moment (t^). As explained in 
3|], the absence of the term 2Ai {t)^^K/Q in ( [iTsD s an artifact of using the single 



gluon approximation. On the other hand the absence of the term \2JS? /Q^ is physically 
meaningful (a SDG analysis is the appropriate tool to identify such a power-correction, 
if it exists). Our current finding that the effect of A2 is cancelled at the level of the 
distribution is consistent with that of eq. ( p.l5| ). It is interesting to note, however, that 
eq. ( |115[ ) was derived in based on the entire characteristic function, while our current 



conclusion concerning A2 is based on the logarithmically enhanced cross-section alone. If 
two-jet dominance applies to (t^), fixing Ai and A2 by the distribution would allow the 
determination of the power-corrections to (t^) based on ( |1 14| ) . Unfortunately, this is too 
optimistic: as one can verify, e.g. by looking at the data for the distribution scaled by 
(see iSlI]), at Mz this observable gets a significant contribution from the region t ^ 0.2. 



As we saw at various stages of the current investigation (see e.g. fig. |I| and table |^) the 
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logarithmically enhanced cross-section ceases to dominate at t ^ 0.2. It follows that the 
second moment of the thrust depends on contributions we neglected here. Perturbative 
and non-perturbative physics of three jets^, including terms that are sub-leading in t at 
NNLO, must be taken into account for a reliable analysis of this quantity. 



5.4 Theoretical uncertainty 

As always, there are unproven assumptions and various limitations to our approach. 
Here we address several issues shortly. A quantitative analysis is restricted to those 
approximations we have made where a particularly significant impact on the extracted 
value of as is expected. 

First in line is the uncertainty in the perturbative calculation. Here there are two 
major issues: the first is the uncertainty due to missing NNLO calculations. The signifi- 
cance of NNLO corrections was emphasized recently in [^. In our analysis it is reflected 



for instance in the large renormalization scale dependence in the region t ~ 1/3. Once 
such a calculation is available, matching with the DGE would allow a significant reduc- 
tion in the uncertainty, in particular, concerning the extraction of a^. The second issue 
concerns exponentiation: the DGE performed here is only approximate beyond NLL. An 
exact NNLL calculation would allow one to test certain aspects of this approach and 
complement the resummation. 

Concentrating on the two-jet region we neglected effects that are explicitly suppressed 
by t, but renormalon factorial growth is not restricted to the logarithmically enhanced 
terms. Large numerical coefficients appear also in the "remainder function" . Fig. |I| shows 
that at t ~ 0.2, such sub-leading terms become important. Eventually, DGE may be 
generalized to included them. Currently it is not known whether these terms exponentiate 
or not (and therefore they are treated differently in different "matching schemes"). One 
can give an estimate of the possible impact of such terms by making a modification of 
the log, as suggested in ||^, e.g. Inl/t — > ln(l/t — 1). Such a modification ensures 
that the logarithmically enhanced cross-section will vanish at t = 1/2, like the physical 
cross-section, rather than at t = 1. Some results for a shift-based fit are summarized in 
table ^ 

The general quality of these fits is worse than in the case of the non-modified log, 
table 0. In addition there is a stronger dependence of as on both cuts. These two facts 
suggest that such a modification is not quite the appropriate way to take into account 
the fact that the physical cross-section vanishes at t = 1/2. Still, as a rough measure of 
the impact of the terms we neglected these results are useful. For the favoured cuts the 
difference in as is quite appreciable. It changes from ~ 0.110 to ~ 0.114. A similar study 
for the shape-function-based fit (with tL(Mz)) shows a somewhat smaller change in as'- 
the central value changes from ~ 0.109 to ~ 0.112. Altogether we see that NNLO (and 
higher order) terms which are not logarithmically enhanced have a relatively significant 
impact on a^. We consider this as the largest source of uncertainty. 

One aspect in which the DGE may be improved is the running coupling formula. Our 
entire analysis was performed with a two-loop running coupling. The comparison with 



''Recently, non-perturbative effects in three-jet observables were addressed 



48 



Table 6: Fits of as and Ai to data with a modified log. The upper limit is tn = 0.32. 
The results are presented as a function of the lower limit. 



tL(Mz) 






^(Mz) 


Ai 


xVdof 


points 


0.02 


0, 


,1100 


± 


0.0003 


2.24 ± 0.10 


1.16 


231 


0.03 


0, 


,1095 


± 


0.0003 


2.42 ± 0.15 


1.08 


213 


0.04 


0. 


,1106 


± 


0.0005 


1.80 ± 0.25 


0.94 


193 


0.05 


0. 


,1119 


± 


0.0006 


1.10 ± 0.26 


0.88 


177 


0.06 


0, 


,1134 


± 


0.0007 


0.42 ± 0.29 


0.83 


165 


0.07 


0, 


,1146 


± 


0.0010 


-0.11 ± 0.36 


0.85 


152 


0.08 


0, 


,1151 


± 


0.0012 


-0.31 ± 0.41 


0.84 


140 


0.09 


0, 


,1154 


± 


0.0014 


-0.43 ± 0.49 


0.87 


127 


0.10 


0, 


,1160 


± 


0.0017 


-0.64 ± 0.58 


0.88 


119 


0.11 


0, 


,1154 


± 


0.0022 


-0.39 ± 0.84 


0.89 


109 



the one-loop case in fig. ^ shows a relatively large variation. We expect the difference 
between three-loop and two-loop to be much smaller, but it may be non-negligible. 

Next, there are uncertainties in the way hadronization corrections are treated. Here we 
assumed that the dominant non-perturbative corrections appear in a way that matches 
the ambiguity in the perturbative calculation. The success this general approach has 
had in event-shapes and elsewhere is encouraging. However, here we take it one step 
further: we assume that the u dependence of each renormalon ambiguity (not only the 
first!) is reflected in the dependence of the corresponding non-perturbative corrections. 
Of course, any assumption we took in the perturbative calculation may have an impact 
on the prediction. In particular, 

• Having neglected effects that are explicitly suppressed by t in the perturbative 
calculation, we cannot expect the power correction analysis to be reliable beyond 
0{l/t) accuracy. 

• Since exponentiation is based on two-jet kinematics, non-perturbative effects that 
are associated with three jets (where the recoil of the quark must be taken into 
account) are not included. 

• The normalization of the thrust variable by the sum of energies (Q) in (|l^) (rather 
than by the sum of the absolute values of the three-momenta) was shown to affect 
the coefficient of the leading power-correction in the case of the average thrust |l| . 
There is a similar impact on the 1/Qt correction in the case of the distribution, 
and probably also on higher power corrections. 

• Having neglected the non-inclusive decay of the gluon in the perturbative calcula- 
tion in sec. 2, the effect is missed also on the non-perturbative level. See discussion 
below eq. (^Sj) . 
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• In our framework it is very natural to treat non-perturbative corrections as if they 
all exponentiate, since the perturbative calculation indicating their necessity is a 
calculation of the exponent. This is most clearly put in eq. (|67|): the ambiguity 
appears in the exponent and so should be the power corrections. It is probably true 
that in the peak region, the leading power corrections indeed exponentiate, but in 
general there are power corrections that do not. 

There are also hadronization effects that cannot be associated with ambiguities in the 
perturbative calculation. For example, the latter does not contain information on the 



spectrum of the hadrons produced. In a recent work Salam and Wicke ||5^ showed that 
hadron mass effects introduce important corrections, for average values of event-shape 
variables. 

Another fact we ignored here is that in reality also heavy quarks are produced. Since 
the data contains all events whereas in the calculation the quarks are treated as massless, 
some systematic error is expected. 



6 Conclusions 

Many infrared safe observables in QCD are sensitive to soft and collinear gluon radi- 
ation. This sensitivity appears at the perturbative level in the form of large Sudakov 
logs. In order to deal with this situation one must perform resummation of the logarith- 
mically enhanced cross-section. The resummation has two aspects: exponentiation and 
integrals over the running coupling. The Dressed Gluon Exponentiation (DGE) differs 
from previous resummation methods in the way the second aspect is addressed. The ap- 
proach makes a direct link between the resummation of Sudakov logs and that of infrared 
renormalons. We showed that a careful treatment of running coupling integrals makes it 
possible to take into account a large class of logs, which is fully consistent with both the 
exact exponent up to NLL accuracy and renormalization group invariance. The latter is 
realized only when all the powers of the log are resummed. Moreover, sub-leading logs 
are factorially enhanced compared to the leading ones. As a consequence the standard 
approximation based on keeping the leading and the next-to-leading logs alone is not 
always justified. 

DGE is a method to calculate a well defined class of log-enhanced terms based on 
exponentiation of the Single Dressed Gluon (SDG) renormalon-sum. The basic idea is 
that the exponent can be represented in terms of an integral over an observable-dependent 
characteristic function times the running coupling, similarly to the standard dispersive 
approach in renormalon calculations |26|, |2^, [^. This representation is natural since 
the exponentiation kernel is primarily associated with a single gluon emission. The 
method is quite general and can readily be applied to other physical quantities. The first 
stage is to calculate analytically the large Nf or massive gluon characteristic function, 
and then to identify the phase-space limits in the integration over the gluon virtuality as 
well as the specific terms in the characteristic function that contribute to logarithmically 
enhanced terms in the perturbative coefficients. This allows to obtain the logarithmically 
enhanced part of the SDG cross-section to all orders. At the next stage the physical cross 
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section is obtained by exponentiating the SDG result. In the case considered here this 
exponentiation is rather straightforward, given the assumption that gluons are emitted 
independently from the quarks and contribute additively to 1 — thrust. In making this 
assumption we ignored correlations in multi-gluon emission that occur at higher orders. 
As a result, the exponentiation is only approximate as far as sub- leading logs (NNLL 
and beyond) are concerned. DGE cannot replace the construction of the exact exponent 
as a solution of the evolution equation up to a given logarithmic accuracy. It is rather 
complementary to it, similarly to the way resummation in general is complementary to 
fixed-order calculations. 

One of the main results of this work is the observation that sub-leading Sudakov 
logs are factorially enhanced compared to the leading logs. This observation was made 
first at the level of the SDG cross-section, before exponentiation. At this level factorial 
growth appears only due to the collinear limit of phase space. The Borel function (in 
the large /3o limit) has simple poles at z = 1 and 2. However, the Borel integral is still 
well defined thanks to the analytic continuation. A much more dramatic factorial growth 
appears at the level of the exponentiated result in the Laplace conjugate variable u. 
Here both the collinear and large-angle (soft) limits of phase space contribute to infrared 
renormalons. As a result of the large-angle contributions, the first infrared renormalon 
appears at z = 1/2, thus enhancing the divergence of the perturbative expansion as 
well as the relative significance of sub-leading logs compared to the leading logs. In the 
standard exponentiation formula P] the most prominent enhancement of sub-leading logs, 
associated with the z = 1/2 renormalon, appear due to the leading term in the splitting 
function, the same term that generates the leading logs. 

The traditional formulation of resummation is based on writing the exponent as an 
expansion in the coupling, with coefficients gniO that are functions of the combination 
^ = LPoas/n, as in (0) (or in (^)). This way of organising the expansion is convenient 
when working with a fixed logarithmic accuracy. In this formulation, the enhancement 
of sub-leading logs translates into the properties of QniO'- these functions become more 
singular at ,^ = 1/2 and ^ = 1 as n increases, and in addition they increase factorially. 
As was shown in fig. ^, the convergence of the expansion ( ^6|) depends strongly on the 
coupling (and thus on Q^). In the case of the thrust sub-leading logs are quite important 
at all relevant energies. When using DGE and keeping all the logs the standard way 
of organizing the expansion loses its attractiveness: as a consequence of the factorial 
divergence of the functions gn{0 ^t large n the series must be truncated. Truncation at 
the minimal term would minimise the error, but this is not practical since the enhanced 
singularity of increasing order functions implies that the truncation order must itself be 
a function of L. An elegant alternative was presented in section. 3.3: the exponent can 
be written in terms of a finite set of analytic functions /^(InQ^/A^) with coefficients fp 
that are polynomials in L. Each function /^(InQ^/A^) corresponds to a single renormalon 
integral over a simple or a double pole. In this way running coupling effects are resummed 
first and the logs only later. As was shown in fig. ^ the convergence of this resummation 
method is quite remarkable. 

The ambiguities of the perturbative result obtained by DGE can be used as a template 
for the non-perturbative power corrections. In particular, we saw that it is possible to 
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identify the t dependence associated with each renormalon singularity in the exponent. 
The leading power-correction for t ^ K/Q scales as 1/Q. It corresponds to a simple 
renormalon pole at 2; = 1/2, with an ambiguity which is proportional to 1/t. This 
results, upon exponentiation, in a shift of the distribution, as predicted by [0, ^ 
The physical origin of this correction is the sensitivity of the thrust to large-angle soft 
gluon emission. The same sensitivity gives rise to other renormalon singularities that are 
further away from the origin. The corresponding ambiguities scale as ~ A"/Q" and are 
therefore less important at large t. However, at small t they behave as l/t*^, implying 
non-perturbative corrections of the form XnK^/{Qt)'^- As emphasised by Korchemsky 
and Sterman, these corrections are all important in the distribution peak region. Since 
they exponentiate, it is natural to describe them by a non-perturbative shape-function. 
The parameters A„ then correspond to the central moments of this function. 

Our calculation shows that the ambiguities ~ A"/(5" with n = 2,4, . . . vanish in the 
large (5q limit. This strongly suggests that the parameters A„ for even n are suppressed 
according to eq. (0). In addition, the appearance of a double renormalon pole at 2; = 1 
leads to a coUinear power correction proportional to A2 A^/ (Qt)^, which has the opposite 
sign to the A2 term in the large-angle contribution. This leads to an almost complete 
cancellation of the overall effect of A2. 

The immediate phenomenological implication of these findings is that by shifting the 
distribution one should be able to fit the data in a wider range (closer to the peak) than 
a priori expected: there is no new effect beyond that of the shift so long as ^ <^ 1, 

instead of ^ ^ -C 1. The phenomenological analysis indeed shows very successful shift- 
based fits, even for rather low t values (this was first observed in 0). Clearly, in the 
peak region itself, a shift is not enough, and a convolution with a shape-function must 
replace it. 

Finally, our main conclusions from the phenomenological analysis can be summarised 
as follows: we obtained good fits based on the DGE result both with a shifted distribution 
and with a shape function. The two methods are fairly consistent with each other. This 
demonstrates that the data in the peak region can be used to extract as provided the 
appropriate resummation is performed and the power corrections are included. In the 
fits, the central value of af^ (Mz) ranges between 0.108 and 0.115. The most important 
factor which limits the determination of is the absence of a NNLO calculation. In 
particular, terms that are explicitly suppressed by t at NNLO are important, as can be 
deduced from the fact that the modification of the log: In 1/t — > ln(l/t — 1) can change 
the extracted ^ (Mz) by 4%. There are also non-negligible differences (3%) between 
different fitting procedures, such as dependence on the low t cut, the use of a shift or a 
shape-function, and the parametric form of the latter. Other factors contributing to the 
uncertainty are discussed in sec. 5.4. 

The correlation between the non-perturbative parameters A„ and as is strong. Con- 
sequently a quantitative discussion is restricted to the case where is fixed. Then the 
shift parameter, or the first moment of the shape function, Ai is reasonably well deter- 
mined (see fig. ^. Even when a modified log is used, the value of Ai as a function of as 
remains roughly the same. The sub-leading power-corrections are harder to constrain by 
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the data. To determine the higher central moments of the shape-function it is useful to 
fix also Ai. As an example, fixing as and Ai to their best fit values we find that A2 is 
small and that the fit is not affected much by fixing it to zero (see table |^). Thus our the- 
oretical prediction that the effect of A2 is small can be consistent with the data. A clear 
preference of the fits for a shape-function which is not positive-definite was observed. 
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A Numerical calculation of the NNLL part of the 
NLO coefficient 



To calculate the NNLL part of the NLO coefficient we have used the EVENT2 |TB 
program. Writing the thrust-distribution on the form 



A{L)a{Q') + B{L)a{Q')\ (116) 



where L = \n{l/t), a = af^ /ir, and A{L) is the analytically known leading order coeffi- 
cient, the program calculates the next-to-leading order coefficient B{L) in the MS scheme 
using Monte Carlo integration. Note that the normalisation in Eq. ( |116D is to the Born 
cross-section o"o and not to the total cross-section a = (l + |Cira(Q^) + . . ao- The 
program is written such that the different colour factor parts can be calculated sepa- 
rately (this was done before using the EVENT program in ^) in the following 
way, 

BiL) = CIBcAL)+CfCaBcAL) + CfTfNfBtAL)- 

Since we are interested in the large L (small t) region we decreased the invariant 
mass cut-off on any pair of partons, mfj > CUTOFFQ^ in the program to 10 from 
the default 10~^ in order not to restrict the phase-space for small t. To increase the 
sampling of the small t region we also changed the parameters NPOWl and NP0W2 
to 4 instead of the default values 2. The calculations presented below are based on 10^ 
events (samplings of the integral). We used a linear binning in L and the bin size was 
chosen to be AL = 0.5. 

For each bin in L the Monte Carlo integration gives the integrated cross- sect ion. From 
these integrals we then subtract the known NLL parts (except the constant contributions 
which will be dealt with separately) and normalise to the bin size giving the following 
integrals for each bin, 

h = f'^^^\B,{L)-Br{L))dL, 
where i = Cp, Ca, Tp, and 




The results for the different integrals are shown in fig. |I2| If it were not for sub-leading 
corrections, the integrals Jj would coincide with the NNLL coefficients Bi^i. In principle 
one could take into account sub-leading terms by fitting a more general form to the 
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Figure 12: Results of the numerical calculation of the NLO coefficients Bi with the 
NLL terms subtracted for the different colour factors. The band shows the estimated 
systematic error from sub-leading corrections when extracting the NNLL coefficients. 
The results of the fits are shown as straight lines in the regions fitted. 



integrals !{. We choose instead to estimate a systematic error from sub-leading corrections 
for each bin in L. The estimate is based on assuming that the most important form of 
sub-leading corrections has the same power of logs as the leading log multiplied by one 
factor of t. The size of the coefficient multiplying this term is then assumed to be the 
same as Bi^i estimated from the integral Jj. The systematic error for each bin is thus 
estimated by, e.g. 

^^c, = -TF / L'exp{-L)dL. 



AL 

These errors are illustrated as bands in fig. |r^. When making the fit these systematic 
errors are added in quadrature to the statistical errors from the Monte Carlo integration. 



From fig. it is clear that for L ^ 7 the sub-leading corrections start to become 



important. It is also clear that for L ^ 12 the statistical errors start to increase signifi- 
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cantlypj. Thus we choose the fit range 7.1 < L < 12.1. Fitting a constant in this range 
gives the following results, 



2.95 ±0.27 ±0.11, 
1.855 ±0.081 ±0.093, 
-1.255 ±0.022 ±0.010. 



where the first error is from the fit and the second from the maximal variation when 
varying the upper and lower limits for the fit by ±1. 

From these results we can now get the NNLL term Q21 in the order contribution 
to In R, 



Qi2L^ + QiiL 



a ± 



^23-^^ ± ^22-^^ ± G21L 



For the Ca and Tp parts these are just the same as the Bi^i but for the Cp part one has 
to subtract the contribution coming from the expansion of the exponent and the different 



normalisation used in Eq. (|116|) , 
the result as 



3)/4. Adding the errors in quadrature and writing 



G21 — Gi3o,iCfPo ± Ocf,iCf + Gca,i^aCf 



we get 



1.23 ±0.29, 
-1.60 ±0.14, 
3.766 ±0.070. 



These results can be compared with Eq. ( |107|) where 



0.13, 

-0.73 

3.75, 



in order to examine to what accuracy the DGE generates the sub-leading logs. The coef- 
ficient which is directly calculated in our approach is Gp^^i- In this case the agreement is 
very good. As expected, the other coefficients which are generated by the exponentiation 
do not agree. It is encouraging, however, that their values are smaller and that they 
have the correct signs. Thus, including the corresponding terms at NNLO and beyond is 
certainly an improvement. 

For easy comparison with earlier results ^ ^ we also write the result in the 
expansion 



Ini? 



Gi2L^ ± GuL 



a/2 + 



G2zL^ ± G22L'^ ± G21L 



aV4± 



^For even larger L the phase-space limitation imposed by the invariant mass cut-off becomes clearly 
noticeable. 
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with G21 = Gtp,i + Gcp,i + Gca,i- Our results (for Np = 5) are then given by 

Gcp,i = 9±2, 
Gca,, = 30 ±2, 
Gtp,i = -16.7 ±0.3. 

The resuh for the sum G21 = 22 ± 3 agrees within errors with earher results, but the 



results for the colour factors Ca and Tp are significantly different The reason for 
these discrepancies can probably be found in the different ranges in L used for the fit, 
namely 7.1 < L < 12.1 and 3.1 < L < 6.7 respectively. 
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